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1 Introduction 



Since the beginning of this century we have had two remarkably successful theories, 
quantum theory and general relativity. Each theory describes its own domain with an 
astonishing accuracy. No experiment has ever been in contradiction with either theory. 
However, history provides several examples that attempts to describe apparently un- 
related phenomena within a common theoretical framework have led to great progress 
in our understanding of nature. A unification of quantum theory and general relativity 
might shed light on fundamental questions concerning the nature and origin of space 
time. 

During the past decades the problem of quantizing gravity has attracted the interest 
of an increasingly number of physicists from different fields. Numerous approaches to 
the problem have been developed. Here we will report of some results of approaching 
the problem keeping to the conventional quantum field theoretical framework. It is well 
known that by following this line of action one immediately comes up against severe 
difficulties. General relativity based on the Einstein Hilbert action is not perturbatively 
renormalizable in dimensions higher than two. Adding higher derivative terms to the 
interaction does not improve the situation. Renormalizability can be obtained but then 
unitarity is lost [I]. In 1979 Weinberg pointed out that under special circumstances one 
can construct a reasonable quantum field theory out of a theory apparently needing an 
infinite number of counter terms [2]. The special circumstances required are that the 
original theory has an ultraviolet fixed point for which the associated critical surface 
is finite dimensional. In that case only a finite number of parameters is needed to 
describe the ultraviolet behaviour of the theory and the situation resembles the one 
encountered for a renormalizable theory. This scenario is known as asymptotic safety. 
Weinberg and others have studied in perturbation theory the renormalization group 
fiow of Newtons constant using dimensional continuation [2]. The outcome of these 
studies is that gravity appears to be asymptotically safe in 2 + e dimensions (e > 0). 
Since the analyses are based on perturbation theory and e is assumed to be infinitesimal 
we can not really take these results as an indicator of what happens in the ultraviolet 
regime of 4-dimensional gravity. Nevertheless Weinbergs ideas are sill viable and the 
work we will report on is greatly infiuenced by them. 

An approach to quantizing field theories which has made possible addressing non 
perturbative questions is the Euclidean path integral method. At the same time as 
Weinberg put forward his ideas. Hawking suggested to apply the path integral ap- 
proach to the problem of quantizing gravity [3]. In particular he suggested that the 
rotation from Lorentzian to Euclidean space in the case of ordinary field theories should 
be generalized to a rotation from Minkowskian to Riemannian space. Writing down 
a Riemannian path integral for quantum gravity is straightforward but the possibility 
of recovering Minkowskian space times from Riemannian ones by some kind of analyt- 
ical continuation remains an assumption. Even though one chooses to postpone the 
clarification of this point, ones troubles are not over. The Riemannian action is un- 
bounded from below due to the conformal model. Several suggestions for a solution of 
this problem exists. Hawking et. al. have suggested that one might rotate the contour 
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of integration for the conformal factor so that it becomes parallel to the imaginary 
axis thereby making the path integral convergent [4]. Greensite has advocated that 
stochastic stabilization provides a way of attributing meaning to Riemannian quantum 
gravity [5]. Here we will analyse the Riemannian path integral from what one might 
call a statistical mechanical point of view. By discretizing the path integral we con- 
struct a lattice model of quantum gravity. We then search the coupling constant space 
of the discretized model for critical points, the hope being that a second order phase 
transition point where physical observables show appropriate scaling behaviour exists. 
If such a situation is indeed encountered we might attempt to define a continuum limit. 
We deal with the problem of the unboundedness of the Riemannian action by choosing 
a discretization which for fixed space time volume and topology provides a regulariza- 
tion of the path integral. Of course this means that in the first place we just postpone 
the problem untill the infinite volume limit is taken. Several other open questions of the 
Riemannian approach to quantum gravity turn up also in the discretized models. One 
example is the question of how to choose the path integral measure. Another example 
is the question of whether or not one should sum over space time manifolds of different 
topology. In two dimensions the sum over topologies is known to be badly divergent 
and the situation only gets worse in higher dimensions, let alone the problem that no 
classification of manifolds according to topology is known in higher dimensions. There 
is of course also still the question if we can extract any information about Minkowskian 
gravity from the Riemannian model. Finally in the lattice approach a new problem 
turns up. General relativity is reparametrization invariant and we want a quantum 
theory of gravity to possess the same symmetry. By discretizing the Einstein Hilbert 
action we break explicitly reparametrization invariance. Therefore it is of outmost im- 
portance in case the lattice theory offers to us the possibility of defining a continuum 
limit to check whether reparametrization invariance is recovered in this limit. 

Despite the non negligible number of unclarified points it is still the hope that some- 
thing interesting can be learnt from studying discretized quantum gravity. Needless to 
say that predictions of experimental relevance are not to be expected in the immediate 
future. It is tempting the repeat the statement made by A. Salam at the first Oxford 
conference on quantum gravity in 1975 [6] and repeated by C. J. Isham at the second 
Oxford conference on quantum gravity in 1980 [7] 

"In particle physics we have become conditioned never to ask what the 
theory can do for us; instead we humbly try to see what we can do for the 
theory." 

In the following sections we will describe some examples of what has been done for 
the theory. In section 2 we introduce the discretization scheme known as dynamical 
triangulations. The remaining sections describe the application of this discretization 
scheme in two, three and four space time dimensions, focusing on aspects addressed in 
the work of the author. 
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2 Dynamical Triangula! ions 



2.1 The Discretization Programme 

The first step in tlie discretization programme consists in replacing tlie continuous d- 
dimensional space time manifold by a special type of lattice known as a J- dimensional 
simplicial manifold. The building blocks of such a lattice are so-called J-simplexes. A 
2-simplex is a triangle, a 3-simplex is a tetrahedron. The J-dimensional analogue can 
be thought of as an object consisting of {d-\-l) vertices connected via (^+^'>^ links. A d- 
dimensional simplicial manifold is a collection of J-simplexes which are glued together 
along their [d— l)-dimensional sub-simp lexes in such a way that the set of J-simplexes 
surrounding any given vertex is homeomorphic to the unit ball, B'^, in R^. Hence 
per definition a simplicial manifold does not have any boundary. In the following we 
will in addition assume that our simplicial manifolds are connected. Let us introduce 
the notation that Ni denotes the total number of z-simplexes in such a J-dimensional 
manifold while rii refers to a given z-simplex. Furthermore let us denote by o[ni) 
the order of n^, i.e. the number of J-simplexes to which rii belongs. Then the Euler 
characteristic of a J-dimensional simplicial manifold, Xd^ is given by 
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Furthermore from the definition of a simplicial manifold it follows that in d dimensions 
the following relations hold between the A^^'s [8] 



For d = 3 and J = 4, (2.2) gives rise to two independent relations between the A^^'s. 
For d = 2 only one constraint arises. In the case d = 3 one of the constraints is that 
the Fuler characteristic has to equal zero. This holds always for simplicial manifolds 
of odd dimension. 



One can introduce the notion of distance on a simplicial manifold, i.e. equip the mani- 
fold with a metric by specifying the lengths of all its links and assuming the metric to 
be fiat in the interior of its d- and [d — l)-dimensional sub-simp lexes and continuous 
across faces. In this way one obtains a so-called piecewise linear space. We note that 
the number of edges of a J-simplex equals the number of independent components of 
the metric tensor in d dimensions. To obtain the discrete analogue of a Riemannian 
manifold restrictions must be imposed on the edge lengths in order to ensure that the 
metric is positive definite. The construction of piecewise linear spaces is due to Regge 
who also showed how one can assign curvature to such spaces [9]. Since the geometry 
of a J-dimensional piecewise linear space is Fuclidean on its d- and [d— l)-dimensional 
sub-simplexes its Riemann tensor is (5-function like having support only on the [d — 2)- 
dimensional sub-simplexes, the hinges. We will not enter into a discussion of how to 
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define the full Riemann tensor for a piecewise linear space. Such a discussion can be 
found in reference [10]. 

Our aim here is only to find the discretized version of the Riemannian Einstein 
Hilbert action given in the continuum by 

Seh[M] = a / d'x^ - / d'x^R (2.4) 

where Ai is the space time manifold, A the cosmological constant and Gjv Newtons 
constant. Hence as far as curvature is concerned all we need to know is how to calculate 
the integrated scalar curvature. The integrated scalar curvature for a given piecewise 
linear space is a sum of contributions from each of its hinges, namely 

d'x^R^ J2 V{h)6h (2.5) 

■M I,- I, 

nznges,n 

where V{h) is the volume of the hinge, /i, and 6^ the deficit angle associated with the 
hinge. The deficit angle at a hinge, /i, is 27r minus the sum over J-simplexes containing 
/i, n^^\h)^ of the angle between the two faces of n^^\h) which have h in common. 
In general the deficit angle 6^ is a complicated function of the edge lengths of all d- 
simplexes meeting at h. A discussion of the general case can be found for instance in 
references [11, 12]. However we will advocate a discretization scheme where one only 
considers piecewise linear spaces built from regular simplexes. For a regular J- simplex 
the angle, O^^ between any two adjacent faces is given by 

cos9d = \ (2.6) 
a 

and hence the integrated scalar curvature for a J-dimensional Riemannian piecewise 
linear space built from regular simplexes can be written as 

/ d'x^R ^ E - oind-2)) = CdNd-2 - Kd+1,2 Nd (2.7) 

where 

Q = ^ and Kd+1,2 =(^'^1^^y (2.8) 

Only for J = 2, is an integer and hence only in two dimensions we can have a 
regular tessellation of fiat space. Furthermore in two dimensions the integrated scalar 
curvature is a topological invariant, namely 

/ d'x^R = i7rx = 87r{l-g), [d = 2) (2.9) 

■JM 

where g is the genus of manifold Ai. The discrete analogue of the cosmological term 
reads 

/ d'^x^^Nd. (2.10) 

JM 



4 



It is also possible to construct discrete versions of integrals involving higher derivative 
terms. However this requires some care since higher derivative terms involve higher 
powers of the Riemann tensor which as mentioned above is (5-function like. We will 
return to this problem in section 4.4.3. From (2.7) and (2.10) it appears that for a 
triangulation, T, consisting of Ni z-simplexes, z = 0, . . . , J, the Einstein Hilbert action 
turns into 

Seh[T] = KdNd-Kd-2Nd-2 (2.11) 

where K,d-2 and k.^ are dimensionless coupling constants, k.^ being a linear combination 
of the cosmological constant and the inverse gravitational constant and K,d-2 being 
proportional to the latter. In two dimensions we can also write 

Seh[T] = K2N2 - /€X, {d = 2). (2.12) 

Having arrived at a lattice version of the action we have completed the first step of the 
discretization procedure. We come now to the second and more subtle step, namely 
the discretization of the path integral measure. In the continuum formulation in order 
to calculate the partition function we are supposed to integrate over all Riemannian 
metrics of our space time manifold modulo diffeomorphisms. This integration we will 
replace by a summation over all simplicial manifolds built from regular simplexes, a 
prescription known as dynamical triangulation ^ 

f d{gM\ Seh[M] , V ^ .-Sp.„\T] /r,-|0\ 

Here we have divided the contribution from a given triangulation by the order of the 
symmetry group of the triangulation, C{T). This can be thought of as eliminating 
the last remainings of the diffeomorphism group. Apart from the symmetry factor we 
choose the same weight for all simplicial manifolds. There are no indications that other 
choices should be made. On the contrary results obtained in two dimensions speak in 
favour of a uniform weight (cf. to section 2.3). One may also argue that some kind 
of universality should hold so that the detailed form of the measure should not be of 
importance. There have, however, been studies of dynamical triangulations with a non 
uniform measure [13]. 

So far we have not made precise what class of manifolds should be included in 
the functional integral. In a theory of quantum gravity one might ultimately want to 
include manifolds of all possible topologies. However at the present stage of studies a 
restriction of topology is necessary. This is evident already in two dimensions. Let us 
write the partition function for dynamically triangulated 2-dimensional gravity as 

Z(«, ^2) = E e'^" |e e"'^^^^^x(iV2) I (2.14) 

X [N2 J 

^ Since we have restricted ourselves to considering only connected simplicial manifolds it is actu- 
ally the free energy rather than the partition function we are calculating by the prescription (2.13). 
However we will carry on the historically conditioned misuse of notation and denote the right hand 
side of (2.13) as the partition function of dynamically triangulated gravity. 
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where 



NAN2)= E T^n- (2-15) 

Ter(N2,x) 

Here J2TeT{N2,x) nie^'iis summation over aU triangulations corresponding to manifolds 
of volume N2 and Euler characteristic x- Two-dimensional dynamically triangulated 
gravity can be solved by analytical means (cf. to section 3.1) and one finds in the limit 
iV2 ^ 00 [14, 15] 

N-^{N2) ~ N]^~^e^'^^', (2.16) 

where 

7^-2 = (7o-2)(l-^/) and 7o = (2-17) 

It is important to note that ^2 does not depend on topology. If we restrict the topology 
the sum (2.14) will be well behaved for k,2 > However, for unrestricted x the sum 
in equation (2.14) will diverge. We can not even hope for Borel summability. In higher 
dimensions the need for a restriction of topology can only be more pronounced. In 
addition we would not even know how to attribute a meaning to the phrase "summing 
over all topologies" since a topological classification of manifolds is not known for 
d > 2. We will take as the partition function for 3- and 4-dimensional dynamically 
triangulated gravity 

^(«.-2,«.) = E -^e-^''^+^^-^''^-^ (2.18) 

where the sum is over simplicial manifolds homeomorphic to S"^ . Considering only 
manifolds homeomorphic to a sphere does not impose further restrictions on the A^^'s 
for J = 3 while for d = 2 and J = 4 it dictates that the Euler characteristic must equal 
two and hence imposes one additional constraint on the A^^'s. In conclusion we see that 
for J- dimensional simplicial manifolds homeomorphic to S'^ we have for d = 2 only one 
independent Ni and for d = 3 and d = A only two independent A^^'s. In particular for 
d = 3 and d = A the action (2.11) is the most general one we can have involving only 
iV.'s. 

Yet no powerful analytical tools for studying dynamical triangulated gravity in three 
and four dimensions have been developed. However, numerical simulations have been 
performed. A priori it is not clear that the number of manifolds with spherical topol- 
ogy and a given volume is exponentially bounded by the volume in higher dimensions 
as it is in two dimensions (cf. to equation (2.16)). However it has now been proven 
numerically both for d = 3 [16, 17] and for d = A [18, 19, 20, 21]. Since for d = 3 and 
J = 4, Nd-2 < const X Nd there exists a line k.^ = n'^^{nd-2) in the '«d) coupling 

constant plane with the property that for Hd > n-din-d-i) the partition function (2.18) 
is well defined. The question of whether the subleading correction to the exponential 
growth is power like in three and four dimensions has also been addressed numerically. 
In three dimensions one finds the subleading correction to be power like only in some 
part of the coupling constant space. The situation might very well be the same in 4 di- 
mensions [22, 23]. We will treat the numerical simulations of quantum gravity in three 
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and four dimensions in more detail in sections 4.3 and 4.4. The dynamical triangulation 
prescription for discretizing the path integral in quantum gravity has many attractive 
features. For a given volume of the space time manifold the Einstein Hilbert action 
is automatically bounded from below as well as from above. Furthermore we have no 
problem with overcounting. Two different triangulations lead to different curvature as- 
signment and hence correspond to two truly different Riemannian metrics. Dynamical 
triangulation as an approach to discretizing the path integral in quantum gravity was 
originally introduced for gravity in two dimensions in references [24, 25, 26, 27]. Let 
us mention that there exists another approach which one could denote as "static tri- 
angulation" which has been advocated by Hamber and Williams. For an introduction 
see for instance [ff]. In this approach one considers a simplicial manifold with fixed 
connectivity and the path integral is calculated by summing over all edge lengths of 
the given triangulation. For static triangulations the path integral is not automatically 
regulated. A cut off on the edge lengths has to be introduced. Furthermore one does 
not automatically have a Riemannian metric. To ensure this further constraints on 
the edge lengths are needed. The question of the choice of the measure is also more 
complicated for static triangulations since many choices of edge lengths correspond to 
the same metric. To divide out the diffeomorphism group is not practically doable. Fi- 
nally one might fear that a fixed triangulation does not allow one to explore completely 
the space of all possible Riemannian metrics. Of course one could also argue that not 
all Riemannian manifolds can be approximated with a simplicial manifold built from 
regular simplexes. However, our task is not to approximate a given Riemannian metric 
but to pick out with appropriate weights representatives from all regions of the space 
of such metrics. Whether we succeed in doing so should of course be tested. Another 
important test concerns reparametrization invariance. As pointed out in the introduc- 
tion, in case the discretized theory offers to us the possibility of defining a continuum 
limit we must make sure that reparametrization invariance is recovered in this limit. 
We will return to these points in section 2.3. Let us just mention that dynamically 
triangulated gravity seems to pass the tests whereas the situation is more doubtful for 
the static triangulations. 

2.2 Coupling to Matter 

A quantum version of general relativity should of course include a description of the 
interaction between space time and matter. Matter fields can easily be coupled to 
dynamically triangulated gravity. The probably simplest type of matter one can in- 
troduce is an Ising spin system. The spin variables, cTj-, can be placed either on the 
vertices or in the center of the J- simplexes of the J- dimensional simplicial manifolds. 
One simply adds the Ising model action 

Sisvag = 1^ (2.19) 
<h3> 

where J2<i,j> denotes the sum over neighbouring pairs of spins, to the Finstein Hilbert 
action and extends the summation over manifolds to include a summation over all 
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possible spin configurations for a given manifold. One can couple any number of such 
spin systems to dynamically triangulated gravity. Other spin systems such as q-state 
Potts models can also be considered. By making use of rules of simplicial differential 
calculus one can couple fields with tensor properties to the geometrical system. For 
example the continuum action for a massless scalar field (f> 

Scout = I d'x^g^''d,<t>d,<t> (2.20) 

translates into 

Slatt^ce = ^ {<Pt " <Pjf ■ (2.21) 
<h3> 

Also in this case we can choose to place the field variables either on the vertices or in 
the center of our J-simplexes. For d = 2 only the model consisting of one Ising spin 
system coupled to dynamically triangulated gravity can be solved analytically [28, 29]. 
Needless to say that for d = 3 and J = 4 no analytical tools for studying dynamically 
triangulated gravity coupled to matter exist. 

Let us try to describe what effects one might expect to encounter when matter is 
coupled to the geometrical system. In case the matter field variables are independent 
as is the case for instance in the disordered phase of the Ising model the geometry can 
not be affected by their presence. However, if the matter field variables are correlated 
they might be able to infiuence the properties of the geometrical system. An example 
of such behaviour is seen in the model consisting of one Ising spin system coupled to 
2-dimensional dynamically triangulated gravity. At the critical point of the spin system 
the critical exponent, 70, describing the fractal structure of space time is changed from 
— I to — |. In dimensions higher than two one could imagine that the presence of 
matter could alter the nature of a possible phase transition in the geometrical system 
(cf. to section (2.3)). The interaction between geometry and matter might also cause 
a change in the critical properties of the matter fields. This effect is seen in the above 
mentioned model as well. Whereas the phase transition of the Ising model on a regular 
2-dimensional lattice is of second order it becomes of third order when the Ising model is 
coupled to two-dimensional dynamically triangulated gravity. In section 4.3.2 and 4.4.2 
we will describe what effects are observed when matter is coupled to dynamically 
triangulated gravity in higher dimensions. 

The prescription given above for the coupling of matter to dynamically triangulated 
gravity is a very direct one. In two dimensions matter can be introduced in a less direct 
way. By considering triangulations consisting of different polygons some of which 
appear with negative weights one can obtain models describing non unitary conformal 
matter with (p^q) = (2,2m — 1) coupled to two-dimensional gravity [30, 31, 32]. These 
models, the Kazakov multi- critical models can be studied by analytical means and are 
the subject of sections 3.2 and 3.3. Recently it has furthermore been shown that by 
coupling an Ising spin system to triangulations of the type just mentioned one can 
construct models describing two-dimensional quantum gravity coupled to any minimal 
conformal model and that by a slight generalization of this idea all (p, q) rational matter 
fields can be reached [33, 34]. Many of the analytical tools developed for the (2, 2m — 1) 
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models are probably applicable also in the generic case. We will comment on this new 
development in section 3.4. The matter fields which can be coupled to gravity in this 
indirect way all have central charge, c < 1, namely 

c=l-^^^. (2.22) 
pq 

For matter fields with c > 1 the direct approach must be used and only numerical 
investigations are possible. Results of such investigations can be found in references [35, 
36, 37]. 



2.3 The Continuum Limit 

The ultimate goal for our studies of dynamically triangulated gravity is of course the 
construction of a continuum theory of quantum gravity. Let us try to describe what 
possibilities for defining a continuum limit the discretized models might offer us. There 
is an important difference between the two-dimensional and the higher dimensional 
cases which can be traced back to the fact that in two dimensions the curvature term 
in the Einstein Hilbert action is purely topological. Let us consider the partition 
function of dynamically triangulated gravity in two dimensions for spherical topology. 
From the expressions (2.f4), (2.f 5) and (2.f 6) it follows that in the limit N2 ^ oo 

Z{n2) \^_, ~ N]°-^e-^^'-^'^^^' (2.23) 
N2 

where we have left out the trivial factor e^*^. We see that for k,2 < ^2 sum will 
be totally dominated by contributions from surfaces with large areas and will not 
converge. On the other hand for k,2 > ^2 surfaces with large areas will be heavily 
suppressed. Only by letting k,2 approach ^2 from above we can arrange that surfaces 
of all sizes contribute to the sum and still the sum is well behaved. Considering the 
limit K,2 — > ('^2)''' hence our only hope for arriving at a continuum theory. In this 
limit we have^ 

^(«2)U~(«2-«D'-"'". (2.24) 

Hence the continuum limit we can obtain in two dimensions is not one to which a 
divergent correlation length is associated. Only moments of the volume of order higher 
than three diverge. In J = 3 and J = 4 we have the possibility of another situation. 
In these cases after having restricted the topology of our space time manifold we still 
have two coupling constants K,d-2 and k.^ at our hand (cf. to equation (2.18)). It follows 
from the discussion in section 2.1 that the partition functions for 3- and 4-dimensional 
dynamically triangulated gravity behave as 

Z{k,_2, «d) = E fi^d, Kd-2)e-^^'-^''^^^'-'^^''' (2.25) 

^Actually also some analytic terms of the type (k2 — H'z)"" , < n < 2 — 70 appear on the right hand 
side of (2.24). This will be evident from the discussion in section 3.2. 
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where f{Nd,Hd-2) is a subleading correction. By the same arguments as used in the 
two-dimensional case one reaches the conclusion that only by approaching the line 
Hd = K-din-d-i) from the region where Hd > n'^^{nd-2) one has the hope of being able to 
address continuum physics. However in this case we still have one free parameter left 
and our hope is that when we move along the line Hd = n-di^^d-i) we will encounter a 
divergence in an appropriately defined susceptibility for quantum gravity, i.e. a second 
order phase transition. This signals the existence of a continuum theory having a 
massless excitation. In practise one does not perform the sum over Nd in (2.18). 
One investigates the behaviour of the partition function for a fixed Nd- However, by 
considering a sequence of increasing Nd^s one can extract information about the infinite 
volume limit. 

In case a second order phase transition point is encountered one must as in ordinary 
statistical mechanics to define a continuum theory prescribe how the critical point 
should be approached. In particular such a prescription must imply a meaningful 
scaling of physical observables. In two dimensions an interesting recipe for defining 
a continuum theory exists [38, 39, 40]. This recipe applies not only to the genus 
zero partition function but to the complete partition function involving all possible 
topologies. It consists in defining a renormalized cosmological constant A^^ by [k,2 — 
K2) = a^^R and sending a to zero and e^'^ to infinity while keeping fixed the parameter 

G-^ = a'e^\ (2.26) 

With this prescription known as the double scaling prescription one can write down a 
power series representation of the all genus continuum partition function where contri- 
butions from surfaces of genus g appear with the factor [G~^ A^j(^Y^~^\ The problem 
of summing over topologies is not solved by this procedure, though, only reformulated. 
The series is still not even Borel summable. However the double scaling idea has re- 
vived the hope that we might eventually be able to address questions in 2-dimensional 
quantum gravity reaching beyond the perturbative expansion. We will return to this 
later in section 3.3 and 3.5. As pointed out in the previous section for d = 3 and J = 4 
we have at the present stage no possibility of addressing the question of summing over 
all topologies. We will describe the present state of the art as regards the possibility of 
defining a continuum theory for space times homeomorphic to the sphere in section 4.3 
and 4.4. 

In case we eventually arrive at some continuum theory we must of course provide 
evidence that it qualifies as a quantum theory of gravity. As mentioned earlier one very 
important requirement to be imposed on such a theory is that of reparametrization in- 
variance. In two dimensions we have the possibility of testing whether reparametriza- 
tion invariance is recovered in the continuum limit. A continuum formalism, Liou- 
ville theory, allows us to study two-dimensional quantum gravity coupled to conformal 
matter fields with central charge c < 1 without breaking reparametrization invariance. 
Liouville theory predicts a set of scaling laws for physical observables, referred to as 
KPZ scaling [41, 42, 43]. One of them is that the contribution to the partition function 
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from surfaces of genus g behaves as 

Zg ~ Ag-"'-'')(i-^). (2.27) 
For unitary matter fields ^str is given by 

Istr = ^ - 1 - V(25-c)(l-c)) (2.28) 

where c is the central charge of the matter field. For non unitary {p^q) matter (2.28) 
is replaced by [32] 

Istr = — r. (2.29) 

p + q - 1 

Liouville theory hence predicts a value for 70 (cf. to (2.24)). For pure gravity, c = 0, 
and the predicted value of 70 is — | . This is exactly what is found using the framework 
of dynamical triangulations. Dynamically triangulated gravity coupled to one Ising 
spin system is the unitary Here the prediction of KPZ scahng is 70 = — 

and as appears from the previous section this prediction is also confirmed. No results 
obtained by dynamical triangulation have ever been in confiict with KPZ scaling. Hence 
it seems that dynamical triangulations do allow us to explore the complete space of 
Riemannian metrics and reparametrization invariance is recovered in the continuum 
limit. Furthermore our choice of the uniform measure is justified. We can of course 
not take the two-dimensional results as a guarantee that reparametrization invariance 
is also recovered in the possible continuum limit of dynamically triangulated gravity 
in higher dimensions. However, since the continuum formalism can not be generalized 
to d = 3 and d = A while dynamically triangulated gravity can, we consider this a 
sound way of continuing the analysis. (We note that dynamical triangulations also 
provide another route for continuing the analysis. While the Liouville approach breaks 
down at c = 1 the coupling of matter fields with c > 1 to gravity can be handled 
in a straightforward way within the framework of dynamical triangulations (cf. to 
section 2.2).) 

Let us close this section by mentioning that the static triangulations have not 
been quite as successful with respect to reproducing the continuum results. In this 
approach coupling an Ising spin system to gravity does not lead to a modification of 
7o as predicted by the formula (2.28) [44]. 
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3 D=2. Matrix Models 



3.1 The Complete Perturbative Solution 
3.1.1 The Hermitian 1-Matrix Model 

It is well known that the partition function of 2-dimensional dynamically triangulated 
gravity given by (2.12) and (2.13) can be identified with the free energy, of the 
following matrix model 

Z[A,iV] = e^[^'^] = / d(Pexp{NTT {--(P^ + X(P^)\ (3.1) 
Jnxn I 2 J 

where the integration is over N xN hermitian matrices and where the gaussian measure 
is assumed to be normalized, if we set 

X = e-^% N = e\ (3.2) 

This can be shown in the following very straightforward manner [45]. The building 
blocks of our two-dimensional simplicial manifolds are equilateral, oriented triangles. 
We can consider the term A Tr ^"^ = A (pij (p^k (pki in the matrix model action as being 
associated with such a building block, t, as indicated in figure 1. 

J 




Figure 1. Labels of a triangle. 

Here we have labeled the vertices of t by z, j and k and equipped its links with matrix 
variables (pi^^ (p^k^ (pki where the order of the matrix indices refiects the orientation of 
the links. The term A Tr (p^ is now the product of matrix elements along the boundary 
of t. With this assignment it is easy to see that one generates by the perturbative ex- 
pansion of F[A, A^] all connected, closed, orientable 2-dimensional simplicial manifolds 
since Wick contractions correspond to gluing together triangles along their links always 
identifying only links with opposite orientation. A given triangulation, T, appears with 
the weight, W{T) 

W{T) = ^L_XN2pfN,-N^+No (3_3^ 

C{T) 

which follows easily from the usual Feynman rules. A comparison of (3.3) with (2.12) 
and (2.13) leads to the identification in (3.2). In the matrix model language the topo- 
logical expansion (2.14) reads [46] 

oo 

F[X,N] = J2N'-'^F,[X] (3.4) 

3 = 
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where the term N'^~'^'^ Fg[\] is the contribution to the free energy from surfaces of genus 
g. Each Fg in (3.4) is a power series in A which as explained in section 2.1 is convergent 
for |A| < Ac = e~'^2. Hence each Fg is an analytic function of A in a region around the 
origin. As also explained in section 2.1 singular behaviour arises when A — )• (Ac)~. In 
this limit one has (cf. to (2.16)) 

F,~(A,-A)(2-'o)M, 70 = "2- (3.5) 

It is the singularity at A = Ac which is interesting from the point of view of continuum 
physics. 

In the following we will consider a completely general hermitian 1-matrix model 
given by 

Z[{5j,iV] = e^[-f^'>'^] = / J(^exp(-iVTr y((^)) (3.6) 



/JVxJV 
where 

w)=f:f<^^- (3.7) 

For the generic model we still have a genus expansion like (3.4). However, the sur- 
faces which are generated by Wick contractions are now built not only from triangles 
but from all types of polygons, and different polygons typically appear with different 
weights. Apart from the partition function we will be interested in expectation values, 
(Q(^)), of operators of the following type 

Q{^) =Ti^^' ... Tr (3.8) 

where we by {Q{cf>)) mean 

{Qm = ^ f d<Pexp{-NTT V{<P)) Q{<P) . (3.9) 

The surfaces which appear by the perturbative expansion of (3.9) have s polygonal 
boundaries of length . . . , A;^. As usual we will only be interested in contributions 
from connected surfaces. In stead of working with the expectation values themselves 
we will work with their generating functionals, the s-loop correlators, defined by 

W{p^,...,Ps) = (tt^-...Tt^--) , s>l (3.10) 

where conn refers to the connected part. With the normalization chosen in (3.10), the 
genus expansion of the correlators reads 

OO -I 

W{p,, . . . ,p.) = E JH^. Wg{p,, . . . ,p,), s>l. (3.11) 

3 = 

In the following we will describe an iterative procedure which allows us to solve the 
model (3.6) completely in perturbation theory, i.e. to calculate Wg{pi^ . . . ^Ps) and 
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Fg for (in principle) any g and any s. In particular we will see how the singular 
behaviour (3.5) emerges in a certain region of the coupling constant space. Likewise 
in other regions of the coupling constant space we will see the emergence of Kazakovs 
m'th multi-critical models characterized by 70 taking the value — ^ [30]. Traditionally 
matrix model are studied by the technique of orthogonal polynomials [15]. However 
as regards explicit calculations away from the double scaling limit this technique is in 
practise only applicable to matrix models with a small number of interaction terms. 
As opposed to this our technique works for a completely general potential. The result 
is expressed without any reference to the original matrix model coupling constants and 
Wg{pi^ . . . ,Ps) as well as Fg depend only on a finite number of parameters, namely 
at most 2 X {?>g — 2 + s) and 2 X {?>g — 2) respectively. The version of the iterative 
procedure that we will present here was developed in [47]. An earlier version can be 
found in [48]. The procedure is based on three ingredients, the loop equations, a loop 
insertion operator, and a set of variables denoted as moments. The loop equations 
and the loop insertion operator are the subjects of the next section while the moment 
description is introduced in section 3.1.3. There exists a similar iterative technique for 
the complex matrix model [49]. Here we will consider only the hermitian case. 

3.1.2 The Loop Equations 

The loop equations are nothing but the Dyson Schwinger equations of the hermitian 
1-matrix model [51, 52]. They consist of a set of relations between expectation values 
which express the invariance of the model under field redefinitions. To derive a version 
of the loop equations appropriate for our iterative scheme it is convenient to consider 
the following transformation of the field, cf) 

^-^+»E,^ = * + ^— . (3.12) 

Under a transformation of this type the measure changes as 

d(f) ^ d(f) (1 + e (tt \ (3.13) 



and we get to first order in e 

2 





iVTr y'((^) e-^Try(^) ^g^^^^ 



Introducing corresponding to the matrix, ^, an eigenvalues density, m(A), with support 
on the real axis [14] this equation can be written as [50] 

- ^W{co) = {W{p)r + ^^(P,P) (3.15) 

where V(cu) = ^2^^! 9j'^'' / j where C is a curve which encloses the support of the 
eigenvalue distribution but not the point u = p. To obtain (3.15) we have also made use 
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of the fact that the expectation value of a product of traces factorizes in the following 

way 

Np 

(Tr <P'^ Tr <P'K . . Tr <P'^) = E IK H Tr <p'")connN'-'''''-' (3.16) 

{P}i=i lePi 

where the sum is over all partitions, {P}, of the set {1, . . . Pi, . . ., PjVp being the 
components of the partition P. It is easy to convince oneself that (3.16) holds by 
considering the effect of gluing together two disconnected surfaces along some polygon. 
In the following we will always assume that u{X) has support only on one interval, 
on the real axis and that 



1 = / u{X)d\. (3.17) 
Jy 

This implies that the 1-loop correlator, W{p), is analytic in the complex plane except 
from a square root branch cut, and that it behaves as ^ as p — )• oo [14]. This is 

the only assumption we need to find the complete perturbative solution of the model. 
Neglecting the last term on the right hand side of (3.15) and making use of the as- 
sumption about the analyticity structure of W{p) we find the following expression for 
H'o(p) 

wM=^i-^\]'-'^^r'T ("8) 

2 Jc 2tti p - Lo - x){lo - y) ] ^ ' 

where x and y are determined by the matrix model potential in the following way 

One strategy for calculating higher genera contributions to W{p) using loop equations 
was outlined by Migdal [52] and later elaborated by David [53]. Equation (3.15) can 
be considered as only the first in an entire chain of loop equations. The n'th equation 
in the chain involves the s-loop correlators for s = l,2,...,n-|-l and can be derived by 
exploiting the invariance of the (n — l)-loop correlator under field redefinitions of the 
type (3.12). As in equation (3.15) the term involving the correlator with the largest 
number of loops is suppressed by a factor Ij N"^ when compared to the rest. This makes 
possible an iterative scheme for calculating Wg{p). However, in each step of iteration 
a new loop equation must be taken into use and calculations quickly become rather 
involved even for simple matrix model potentials. By means of this method W\{p) was 
calculated for even potentials of degree four and six in reference [54]. 

The iterative scheme referred to in the previous section only involves the first in the 
chain of loop equations, i.e. equation (3.15). We avoid the use of other loop equations 
by introducing a so-called loop insertion operator, -j^^y defined by 



d _ ^ ] d , 
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As indicated by its name this operator allows us to step from the n-loop correlator to 
the (n + l)-loop correlator. More generally we have 

d d dF 

At first sight it seems as if nothing is gained by introducing an operator having an 
infinite number of terms. However, as will appear from the next sections, in the 
moment description the loop insertion operator in all cases of interest has effectively 
only a finite number of terms. We note that in order to have a loop insertion operator 
it is necessary to work with a generic potential. A symmetrical potential will not do. 
Now writing W{p,p) as W{p) and inserting the genus expansion (3.11) in (3.15) 
we get 

k - 2Wo{p)} W,{p) = J2 WAp) W,.Ap) + ^W/,_i(p), g > 1 (3.23) 

where K is a linear operator, namely 

A7(p)-/^'^/M. (3.24) 

JC ZTTl p — UJ 

Equation (3.23) expresses Wg{p) in terms of Wg^{p) with gi < g. This equation allows 
us to determine Wg{p) for any value of g starting from the expression (3.18) for Wo{p). 
From W{p) all correlators involving more loops can be found by application of the loop 
insertion operator. Furthermore we can find the free energy, F, from W{p) by, in a 
sense to be explained later, applying the inverse loop insertion operator. 

3.1.3 The Moment Description 

To characterize the matrix model potential it proves convenient to introduce in stead 
of the coupling constants, {(/j}, a new set of variables, {J^, Afj}, denoted as moments 
and defined by 

M,{x,y,{g,]) = i^.- — , k>-l, (3.25) 

Jc 27rz [u — x)'-''^ [u — 

We note that the moments and Mk depend only explicitly on coupling constants, 
(/j, with j > k -\- 1. The moments facilitate the study of the 1-matrix model in several 
ways. First for each term in the genus expansion of the correlators and the free energy 
the dependence of the entire set of coupling constants, {(/i}, arranges into a simple 
function of a finte number of the moments. (Cf. to section 3.1.4.) Furthermore the 
moments refiect more directly than the couplings the possible critical behaviour of the 
matrix model. This can be seen from the following rewriting of Wo{p) 

Woip) = ^V'ip) -^{p- xfl\p - yfl' Y.{iP- ^r'M, + {p- yy-'J,} . (3.27) 

9 = 1 
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To derive equation (3.27) one first deforms tlie tlie contour of integration in (3.18) 
into two, one wliicli encloses the point u = p and one which encircles infinity. In the 
contribution from the latter on then writes [p — as 

1 = ^ [ ^ ] + 1 [ 1 1 (3.28) 

p-u 2\{p- x) - {u - x)] 2\{p-y)-{uj -y)] 



and expand in powers of [J^jz^j and [jj^j respectively. The expansion procedure is 
justified by the fact that the contour of integration encircles infinity which allows us to 
assume that u » p. From the classical work of Brezin et. al. [14] it follows that with 
Wo{p) given as in (3.27) the eigenvalue distribution of the matrix model on spherical 
level reads 

1 oo I 

«o(A) = ^ E {(A - ^y~'M, + (A - yr'J,} ^{x - A)(A -y), y < \ < x. 

(3.29) 

Usually this function vanishes as a square root at both ends of its support. Critical 
behaviour arises when additional roots of Mo(A) approach x oi y. The type of critical 
behaviour is determined by the number, n, of extra zeros at the endpoint where most 
zeros accumulate. Kazakovs m'th multi-critical behaviour is found in the vicinity of 
points in the coupling constant space where n = m — 1. The simplest way of reaching 
a m'th multi- critical point is by demanding that (m — 1) extra zeros accumulate at 
one endpoint, say x, while no extra zeros accumulate at y. The condition for being at 
such a point reads 

Mo = Ml = . . . = M„_i = 0, M„ ^ 0, Ji 7^ (3.30) 

which can easily be read off from (3.29) bearing in mind that the two bracketed terms 
are actually identical (cf. to equation (3.28)). The simplicity of (3.30) is an appeal- 
ing feature of the moment description. An even more appealing feature is that the 
moments as opposed to the coupling constants have definite scaling properties in the 
vicinity of the m'th multi-critical points. This fact makes the moments well suited for 
addressing the continuum limit. We will analyse the continuum limit in sections 3.2 
and 3.3. We note that the expression (3.27) allows us to determine very easily the in- 
verse transformations gi = gi{x, j/, {Mi], {Ji]) since we know that Wo{p) only contains 
negative powers of p. In the moment description the loop insertion operator reads 

d ^ dMn d dJj d dx d dy d 

dV{p) " ^ dy{p) dMn ^ Y dVip) dJj ^ dVip) dx ^ dVip) dy ^' ^ 

where 
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+(n + l/2)A/..„-^ (3.32) 

^ = ^(P-)-'^(p-.)-'^ (3.33) 

Of course the expressions for dJnj dV{p) and dy/dV{p) appear from (3.32) and (3.33) 

by the replacements x ^ J ^ M. To obtain dMn/ dV{p) and dx/dV{p) we have 
rewritten the loop insertion operator in the following way 

d d dx d dy d , , 

+ ^TTTT^ + ^TT^T^ (3.34) 



dV{p) dV{p) dV{p) dx dV{p) dy 

where 



dV{p) fr{ p^+'^ dg 
and made use of the following relation 



V'iui) = 4- — (3.36) 



dV[p) dp p — to 

as well as the boundary conditions (3.19) and (3.20). 

As mentioned earlier to have a loop insertion operator it is necessary to work with a 
generic potential. However, the expressions obtained for correlators and the free energy 
away from the scaling limit can easily be translated to the symmetric case. One simply 
sets X = —y = y/z in the final expressions. This implies setting 

h = (-l)'+'Mfc. (3.37) 

Of course results obtained under the assumption of being close to a critical point of the 
type given by (3.30) can not be taken over to the symmetric case since the possibility 
of having different behaviour at the two endpoints of the cut exists only in the generic 
case. For a symmetrical potential the m'th multi-critical points are characterized 
by the eigenvalue distribution acquiring (m — 1) additional zeros at both ends of its 
support and the condition for being at such a point is the vanishing of the first (m — 1) 
moments. In the vicinity of these points one effectively has a loop insertion operator 
in the symmetrical case too. We refer to [47] for a discussion of this point. 

3.1.4 The Structure of the Solution 

Below we will present the structure of Fg and Wg{p) for the generic hermitian 1-matrix 
model in the moment description. That the structure stated is actually correct to all 
orders in the genus expansion can be proven by induction by means of our iterative 
scheme. We will not go through the proof here. However the course of the proof will 
be evident from section 3.1.5 where we describe the iteration process. 

The genus g contribution to the free energy of the generic hermitian 1-matrix model 
takes the following form 

Fg= j:{a,...a,J,... A /3, 7). " " ,f " ' ^ 9>2 (3.38) 
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where d = x — y is the distance between the endpoints of the support of the eigenvalue 
distribution, ( ) denotes a rational number and a, /3 and 7 are non-negative integers. 
The indices cti, . . . , a^, . . . , /3; take values in the interval [2, 3(/— 2] and the summation 
is over sets of indices. In particular Fg depends on at most 2 X [3g — 2) moments. 
Furthermore, the following relation holds 

F : {ai . . .as; /3i . . . /3i\a,/3,-f) = {-iy'{/3i . . . /3i;ai . . .as\/3,a,-f) . (3.39) 

This follows from the fact that Fg must be invariant under the interchange of x and y 
since nothing allows us to distinguish between the two. In addition there are certain 
restrictions on the integers which enter equation(3.38), namely 

s-a < 0, /-/3<0 (3.40) 

and 

Fg-. {s-a) + il- l3) = 2-2g, (3.41) 

Fg-. ^(a,-l) + ^(/3,.-l)+7 = 45-4. (3.42) 

8=1 J=l 

The relation (3.41) follows from the fact that the partition function Z = e^a " 
is invariant under simultaneous rescalings of N and the eigenvalue density, m(A); 
N ^ k ■ N , u{X) — > |m(A). The relation (3.42) follows from the invariance of Z un- 
der rescalings of the type N ^ ■ — )• y~^gj. Finally the following inequality is 
fulfilled: ^ 

Fg-. ^(a,-l) + ^(/3,.-l)<35-3. (3.43) 

8=1 J=l 

In combination with Fq. (3.42) it gives 

5-l<7<45-4. (3.44) 

The quantities (cti . . . a^; /3i . . . /3, 7)^ can be given a geometrical interpretation. 
We will return to this point in section 3.3. 

To explain the structure of Wg{p), let us introduce the basis vectors x^"'H.p) 
$W(p): 

/''Hp) = ^|$i'^)(p)-Ex('Hp)i^-fc+i}, n>l, (3.45) 
¥-^\p) = jl<^^^\p)-J2'^^'Hp)Jn-k+i], n>l (3.46) 



k=i 



where 



^i'^Hp) = (p-x)-^ {{p-x){p-y)}-'/\ n>0, (3.47) 
^i'Hp) = (p-y)"' {{p-x){p-y)}-'^\ n>0. (3.48) 
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It is easy to show for the operator K defined by equation (3.24) that 
{k-2W,{p)]x^-\p) 

and that the kernel of — 2Wo(p)} is spanned by = 

In accordance with relation (3.22) and (3.38) the genus g contribution to the 1-loop 
correlator can be written as 

3a-i 

WAp) = E {^^Ah^''\p) + D^;)^>^-){p)] , g>l. (3.51) 

n=l 

We note that $(,°)(p) = $^°^(p) does not appear in (3.51). The presence of such terms 
would contradict the boundary condition W{p) — > - for p — )• oo since this behaviour 
was already obtained for genus zero. Furthermore we note that the structure (3.51) of 
Wg{p) is in agreement with the assumption that W{p) is analytic in the complex plane 
except for a square root branch cut [y^x] on the real axis. 

The coefficients A^") have the same structure as Fg (cf. to equation (3.38)) and the 
relation (3.40) is valid also for A^^h However the indices cti, . . . , a^, . . . ^ fii now take 
values in the interval [2,3(/ — n\. The invariance of the partition function under the 
rescalings described above has the following implications for the structure of A^"); 

: {s-a) + {l-fi) = 2-2g , (3.52) 
: j^ia. - 1) + E(/3, - I) + = ^g - 2 - n . (3.53) 

8=1 J=l 

We also have an analogue of (3.43) for A^^\ It reads 

4"^^ ■■ - 1) + - 1) < 35 - n - 1 . (3.54) 

Combining (3.53) and (3.54) one gets again the bound (3.44) on 7. The coefficient 
_D(") appears from A^") by the replacements d — )• —J, J M. This follows from the 
fact that Wg{p) must be invariant under the interchange of x and y. (We note that 
we do not have a relation like (3.39) for the A^^^'s.) Summing up the information just 
given about Wg{p) one finds that Wg{p) depends on at most 2 X [3g — 1) moments. 
Then taking a look at the loop insertion operator (3.31)-(3.33) bearing in mind the 
relation (3.22) one easily convinces oneself that Wg{pi^ . . . ^Ps) depends on at most 
2 X [3g — 2 + s) parameters. 

There are two cases in which the correlators and the free energy depend on less than 
the maximum number of parameters. One case is when one considers a matrix model 
potential with a finite number of terms, i. e. a potential for which gi = for i > k. Then 
Mj and Jj vanish for j ^ k — 1. The other case is when one considers a symmetrical 
potential. Then the two sets of moments become identical and Wg{pi^ . . . ,Ps) depends 
on at most {?>g — 2 + s) parameters, the free energy on at most {?>g — 2) parameters. 



1 

[p — xy 
1 

{p - yY 



n > 1 
n > 1 



(3.49) 
(3.50) 



20 



3.1.5 The Iterative Procedure 



According to Eq. (3.23), in order to start the iterative procedure we need to calculate 

d 

dV(p2) 



Wo{p,p). To do this we first find Wo(pi,P2) by applying ,y'i ^ in the form (3.34) to 



Wo{p) written as in equation (3.18). This gives [55] 

^ 1 I P1P2 - lipi + P2){x + y) + xy 1 

Wo{pi,p2) = —, riT < I = - 1 ? • (3.55) 

2(pi - P2? { ^{p, - x){p^ - y){p2 - x){p2 - y) J 

Then taking the limit pi — )• p2 we obtain 

Now we can determine W\{p) and we see that it is of the form (3.51) with 

D[- ^ 1, . 1. ,3.58, 

Carrying on the iteration process is straightforward. In each step one must calculate 
the right hand side of the loop equation (3.23). Decomposing the result obtained into 
fractions of the type [p — x)"", [p — y)~^ allows one to identify immediately the coeffi- 
cients and (cf. to equations (3.49), (3.50) and (3.51)). To calculate Wg{p,p) 
one makes use of the expression (3.31)- (3.33) for the loop insertion operator. We note 
that there is no simplification of the algorithm in the case of a symmetric potential. We 
can only put x = —y = ^/z at the end of the calculation. This complication of course 
stems from the fact that we have to keep the coupling constants with odd indices in the 
loop insertion operator until all differentiations have been performed. Only hereafter 
can they be put equal to zero. The same is not true in the double scaling limit however. 
This point is explained in reference [47]. 

By taking a closer look at the loop insertion operator (3.31)-(3.33) and bearing in 
mind the results (3.56), (3.57) and (3.58), it is easy to convince oneself that A^") and 
_D(") depend only on x and y via (x — y) and have the structure shown in equation (3.38). 
The results for g = 2 obtained with the aid of Matbematica read 

^(1) ^ 201 67 J2 5J3 49 

^ 256 J5 J^2 ^28 J4 J^3 32 j^s + 256 

57 IIJ2 49M22 



64 J5 Ml 128 J4 J^2 256 M^^ 

201 22 M2 J2M2 



256J5Mi2 128 d^JiMi"^ Gid^ J^"^ Mi' 

67 M2 5M3 
128 J4 Ml 3 ~ 32d3M^ ' 
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57 8J2 49 

201 3M2 J2M2 



+ 



+ 



256 J4 Ml 2 
67 M2 
"l28 JSMi^ 



256 J Mi^ 32JMi3 128J2Mi 



128 J3 Ji Mi^ 

32J3Mi3 ' 
5M3 67 M2 



128 J2 J^"^ Ml' 



201 



+ 



15 



256 J3 Ml 2 128J3JiMi 



5J2 
128 J2 J^^Mi 





49 M2 


189 




128 JMi^ 


256 J2 Ml 


2 ' 




105 






256 JMi^ ' 








= A^') (M ^ 


J, d — > - 


-d) 




= A^') (M ^ 


J, d — > - 


-d) 




= Af^ (M ^ 


J, d — > - 


-d) 




= A^ (M ^ 


J, d — > - 


-d) 




= A^') (M ^ 


J, d — > - 


-d) 



The genus two contribution to W{p) is now determined by equation (3.51). 

Determining Fg when Wg{p) is known is straightforward. The strategy consists in 
writing the basis vectors x^"'H.p) ^^"^(p) as derivatives with respect to V{p). It is 
easy to verify that the foUowing relations hold 



¥^\p) 



dx 



dV{p) 

dy 
dV{p) 

d 

dV{p) 

d 

dV{p) 



2 1 
-In Ml In J 

3 3 

■ - In Ji In J 

3 3 



(3.59) 
(3.60) 
(3.61) 
(3.62) 



Combining this with the results (3.57) and (3.58) one immediately finds 



1 



In Ml 



1 



In Ji 



1 



In d. 



(3.63) 



24 24 6 

For g > I things are not quite as simple. The basis vectors can not be written 
as total derivatives. This is of course in accordance with the fact that the A and D 
coefficients now have a more complicated dependence on the potential. However, a 
rewriting of the basis vectors allows one to identify relatively simply Wg{p) as a total 
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derivative. In the case of x^"'H.p) this rewriting read 



s 




dV{p) i 



J2 X^'^M^.k+i , n>2 (3.64) 

k=2 ) 



where should be written as 



2n-l^^ ^ I " dV{p)i 




.(^) ^ M.^. (3.66) 

The basis vector x^^Hp) should of course still be written as in (3.59). The rewriting 
of the vl/(")'s is analogous to that of the x^^^'s. It can be obtained by performing the 
replacements J ^ M and x ^ y in the formulas above. 

By means of these rewritings we have been able to determine ^2- The result reads 

^ _ 119 119 I8IJ2 181 M2 

3J2 3M2 11J2^ 11 M2^ 



43 M3 43 J3 J2M2 17 

^ 192 Mi^d^ ^ 192 Ji^ J2 + elj^Mi^P ~ 128J1M1 J4 

21J2^ 29J2J3 35 J4 21 M2^ 



160 Ji^ J 128 J/ J 384 Ji^ J 160 Mi^ J 
29M2M3 35 M4 
^128M/J 384Mi3j- ^"^'^'^ 

It is obvious from the formulas above that Fg depends for a non-symmetric potential 
on at most 2 X [3g — 2) moments, and for a symmetrical potential on at most [3g — 2) 
moments. 

3.2 The Double Scaling Limit 
3.2.1 The Partition Function 

The regions in the coupling constant space which are interesting from the point of 
view of continuum physics are those where the free energy becomes singular. From 
equation (3.38) this is seen to take place when Afi(a;, j/, {gi}) or Ji(a;, j/, {gi}) acquires 
a zero of some order. If Mi has a zero in x of order (m— 1) for (x, j/, {gi}) = (xc, j/c, {di})} 
Mi for i < m will at this point have a zero in x of order (m — i). This corresponds 
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to the situation where (m — 1) additional zeros have accumulated at the endpoint 
X of the eigenvalue distribution (cf. to section 3.1.3). If in addition we assume that 
Ji(a;c, j/c, {di}) 7^ we have exactly the situation described in equation (3.30). Let us 
try to analyse the behaviour of Fg in the vicinity of such a singularity. If we perturb 
the coupling constants, {(/i}, away from their critical values, {(/j^}, by an amount, 
Sgi ~ (9(a™), in general [j/c, 3;c] will move to If we assume that 

X — Xc ^ a (3.68) 

we will have 

Mfc~a"-^ ke[0,m] (3.69) 
while the J-moments do not scale. From the boundary condition (3.19) it follows that 

J/-J/c~a". (3.70) 

(An explanation of this point will be given shortly.) A given term in the expan- 
sion (3.38) of Fg [g > 2) under these circumstances scales with a negative power of a, 
P- , given by 

s s 

P_ = a{m - 1) - J2{m - a,) = [a - s){m - 1) + J2{a, - 1). (3.71) 

8=1 8=1 

From equation (3.40), (3.41) and (3.43) it follows that the terms which dominate in 
the limit a — )• are terms for which 

s 

(^a-s)=2g-2, J^^a, - 1) = 3g - 3. (3.72) 

8 = 1 

These dominant terms hence do not depend on any J-moments and have 

1 = 9-1. (3.73) 

In particular we find that 

p-r''{Fg) = {g-l){2m + l). (3.74) 

In deriving (3.71) we implicitly assumed scaling for all Af-moments involved. By a 
slight modification of the argument given above it is easy to convince oneself that 
terms which contain a moment Mk with k > m will be subdominant in the limit 
a — )• 0. Hence in this limit Fg can be written 



as 



3 



F(-)= E ("-•-^H^'lJ^' (3.75) 



l<i:x^<m 



^Here and in the following we use the notation that a superscript (N S) refers to a non symmetrical 
model where the critical behaviour is associated with the endpoint x of the eigenvalue distribution 
while a superscript (S) refers to a symmetrical model. 
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where the indices fulfill the requirements (3.72). Each term in (3.75) is (9(a(^~-^)(^™"'"^)). 
We now see the possibility of the double scaling limit emerging (cf. to section (2.3)). 
Namely, if we send a to zero and N to infinity keeping fixed the product 

G-i = a^^+^N^ (3.76) 

we can rewrite the contributions from surfaces of genus > 2 to the free energy as a 
power series in G. Genus zero and genus one are special. It appears that in the double 
scaling limit Fg for g > I depends on at most [3g — 2) parameters. 

It is easy to convince oneself that if, in addition to the m zeros at the endpoint x, n 
zeros accumulate at the endpoint j/, the scaling behaviour of Fg is still given by (3.72)- 
(3.75), as long as n < m. The case n = m is special, however. In this case Fg is a sum of 
two terms of the type (3.75), one involving Af-moments and one involving J-moments. 
All terms which mix M- and J-moments are subdominant in the double scaling limit. 
This is easily seen to be true for g = 2 from equation (3.67) and is proven in the general 
case in reference [47]. If the potential is not symmetric the double scaling limit of Fg 
for n = m will depend on 2 X [3g — 2) parameters. If the potential is symmetric the 
two sets of moments are identical and only [3g — 2) parameters appear. 

What is traditionally referred to as m'th multi-critical behaviour is obtained when 
the critical point given by (3.30) is approached along the straight line gi = g ■ g^ in the 
coupling constant space. At a given point on this line, assuming the support of the 
eigenvalue distribution to be [y^x] the boundary conditions (3.19) and (3.20) reduce to 

c^{x-x,y^gM^^ + c^{y-y,)gJ^ = 0, (3.77) 
c^{x-x,y^x,gM^^ + c^{y-y,)y,gJ^ = 2{1 - g) (3.78) 

where we have kept only leading order terms and where 

= M,(x„ j/„ {gn) J- = Mxc, yc, {gn) (3.79) 

and 

. = (3.S0) 

From equation (3.77) we get 

{y - vc) = -^-^^{x - xr (3.81) 

■-'1 

which we note explains the statement (3.70). Inserting (3.81) in (3.78) one finds 

{l-g) = ^^c^d,gM^^{x-x,y^ (3.82) 

and with the scaling (3.68) of x we are led to introduce a renormalized cosmo logical 
constant A^^ by 

(l-5)=a"A«. (3.83) 
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Then we obviously have that the Af-moments scale as 



Mfc ~ (aA)/"r-'=M4, A;G[l,m] (3.84) 
and we see from (3.75) the emergence of the formula (2.27) with 

Istr = - — • (3.85) 
m 

It is when the m'th multi-critical point is approached as just described that the specific 
heat (as a function of the renormalized cosmological constant) obeys the celebrated 
string equation [39, 40, 56, 57]. One can consider a more general way of approaching 
the critical point, however. In stead of having only one scaling variable, A^^, one can 
introduce m different scaling variables, /ii,...,//^ corresponding to the m different 
moments of relevance — for example in the following way 

Mfc = a"-Vfc,^ e (3.86) 

With this prescription one finds 

The idea of approaching the m'th multi-critical point in a more general way will allow 
us in section 3.3 to establish in a very direct way the connection between the double 
scaling limit of the partition function of the hermitian I-matrix model and the r- 
function of the kdV hierarchy. 



3.2.2 The Correlators 

In order to be able to speak about the double scaling limit in the context of correlators 
we must introduce a scaling variable tt to replace p. In the case of a m'th multi-critical 
model where the critical behaviour of the matrix model is determined by the conditions 
at only one endpoint of the eigenvalue distribution, x, we set 

p, = Xc + a-K,, (3.88) 
X = Xc + aA^"". (3.89) 

The scaling appropriate for y depends on the conditions at the other end of the eigen- 
value distribution as well as of the details of the prescription for approaching the critical 
point. It must be determined from the boundary equations in each specific case. How- 
ever it always holds that [y — yc) = o[x — Xc). By inserting the scaling relations for 
p, X, y and Af,- in the general expression (3.51) for Wg{p) it is possible by means of 
the equations (3.40) and (3.52)-(3.54) to determine P™"^(Vl/^p(p)) and to derive a set 
of requirements like (3.72) that a given term of Wg{p) must fulfill in order not to be 
subdominant in the limit a — )• 0. Likewise one can study the behaviour of the multi- 
loop correlators by inserting the various scaling relations in the explicit expressions for 
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these correlators. However, by analyzing in some detail the loop insertion operator a 
more direct way of studying the properties of the dominant contributions to the cor- 
relators can be found. From the point of view of the limit a — )• the effect of a given 
operator in d/dV{p) when acting on an expression like (3.38) is to lower the power of 
a by some amount. Examining carefully each term in d/dV{p) shows that the power 
of a is maximally lowered by m + 3/2. All operators which do not lower the power 
of a by this amount hence only give rise to terms which are subdominant in the limit 
a — )• 0. One finds that in the this limit the loop insertion operator reduces to 

d dMn d dx d 

dVip), " ? dVip) dMn ^ dVip)dx ^'^■^■^■^ ^^'^^^ 

where 

^ = -in + l/2)[¥^^'\p)-^mP)} id.s.L). (3.9f) 

/7 T 1 

<5lip) id.s.l.) (3.92) 



dV{p) Ml 



and 



= _ {4 (p - x)}-'/' {d.s.L). (3.93) 



By applying this version of the loop insertion operator to the double scaling relevant 
part of Fg^ i.e. the expression (3.75) (with dc replaced by d) we get the directly the 
dominant contributions to the multi-loop correlators. It appears that in the limit a — )• 
the deviation of y from j/c will never be of importance, so we actually do not need to 
know the scaling behaviour of y. Furthermore we see that all dependence of the J- 
moments disappears in the limit a — )• and that Wg{pi^ . . . ,Ps), </ ^ 1 depends on at 
most {?>g — 2-\- s) moments. From the expression (3.55) it follows that this statement is 
also true for = when s > 2. We note that in the derivation of the expression above 
for the loop insertion operator we implicitly assumed scaling for all moments involved. 
As in the case of Fg it holds for the correlators that for a m'th critical model any term 
which contain Mk with k > m will be subdominant in the limit a — )• and can be left 
out. From the scaling behaviour of the loop insertion operator it follows that 

p^--{Wg{pr,...,p,)) = {m + l/2){2g-2 + s)+s, g > 1. (3.94) 

This shows that for the contributions to the s-loop correlators from surfaces of g > 1 
we need a multiplicative renormalization by a* in order to have a double scaling ex- 
pansion (cf. to equation (3.1 1)). For g = additional renormalization is needed. 
From (3.18) we see that Wo{p) has a part that does not scale as in equation (3.94), 
namely ^V'{p). This part must hence be subtracted. So far we have deliberately not 
used the word "continuum correlators". The precise definition of this concept will be 
given in section 3.3. As we shall see, when we want to make contact with contin- 
uum physics we must also subtract the term — | ^^^^''^^p from Wo{pi,P2) even though 
it has the right scaling behaviour. This is not surprising, though, since — | J^^^p is 
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exactly ^y^^^) \ ^'{Pi)^ the remnant of the term which is subtracted from Wo{pi). 
No additional renormalization is needed in the case of Wo(pi, ■ ■ ■ tPs) for s > 3. 

The scaling prescription (3.88) must be modified if at the critical point the same 
number of zeros accumulate at the two endpoints, x and j/, of the eigenvalue distribu- 
tion. In the case of a symmetrical potential we have x = —y = ^/z and it follows from 
the assumption about the analyticity structure of Wg{j)) that all correlators depend 
only on p via p^. For a symmetrical potential we set 

•p\ = + 2a7riA/5^, z = + 2aA)/™A/5^. (3.95) 

(The normalization of (3.95) is chosen in order to facilitate comparison with the non 
symmetric case.) As mentioned earlier, in the double scaling limit one effectively has a 
loop insertion operator in the case of the symmetrical potential also. We refer to [47] 
for the explanation. Here we just note that as in the asymmetric case Wg{p\, . . . ,Ps) 
depends on at most {?>g — 2 + s) moments except when = and s = I. Similarly, 
what was stated above about the scaling behaviour of the correlators and the need 
for renormalization applies to the symmetrical well. In the very special case 

where at the critical point the same number of zeros accumulate at x and y and yet 
the potential is not symmetric, it is not possible to devise any sensible scaling for p. 
In the following we will refrain from considering this situation. 

3.2.3 The Perturbative Solution 

Having calculated Wg{p) and Fg using the iterative procedure described in section 3.1.5 
one can easily by means of the scaling relations for x, j/, p and the moments determine 
which terms survive in the double scaling limit. However, this is a rather uneconomic 
method for obtaining the double scaling relevant terms. For example from the long list 
of complicated expressions for the A and D coefficients of W2{p) (cf. to section 3.1.5) 
only A^2 \ the first term of ^2^^ and the two first terms of ^2"^^ survive in the double 
scaling limit when we consider a situation described by (3.88) and (3.89). In refer- 
ence [47] an algorithm which gives as output only double scaling relevant terms was 
described. This algorithm not only saves us the trouble of the scaling analysis but also 
allows us to probe higher genera. To develop the algorithm we carried out a simple but 
careful scaling analysis of the loop equation. We do not intend to repeat this analysis 
here. Let us just give a brief outline of the course of the iteration process. 

Let us consider first the case where the critical behaviour of the matrix model is 
determined by the condition at only one endpoint, x, of the eigenvalue distribution. It 
appears from the expressions (3.90)-(3.93) and (3.75) for the double scaling relevant 
part of the loop insertion operator and the free energy respectively that in the double 
scaling relevant part of Wg{p) only the basis vectors x^^\p) will appear and only in a 
form where [p — y) has been replaced by dc- Therefore in the iteration process we must 
use the basis vectors x^^\p) iii this slightly modified form. Apart from this modification 
the recipe is simple. We do exactly as described in section (3.1.5) just starting from 
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the double scaled version of Wo{p^p) 



^i^'Hp,p) = ^(^-^ id.s.L) (3.96) 

and using the double scaled version of the loop insertion operator. Following this recipe 
we have calculated Wg{p) as it looks in the double scaling limit for = 2, 3 and 4. The 
results can be found in reference [47]. 

The procedure for calculating Fg when Wg{p) is known given in section (3.1.5) can 
easily be adjusted to the double scaling limit as well. The starting point is of course 
the double scaling relevant version of Wg{p) — and the strategy consists as before in 
rewriting the basis vectors in a form which allows one to identify Wg{p) as a total 
derivative. However, this time the rewriting in the case of x^") is made with the aid 
of (3.91) instead of (3.32). For genus one we immediately find 

Fi^^^) = -^InMi {d.s.L). (3.97) 

The outcome of the iteration process for = 2, 3 and 4 is 

p(NS) ^ -21 29M2M3 35 M4 
^ 1604Mi^ 1284Mi^ 3844Mi3' ^' ' 



2205 Ma"^ 8685 Ma^Mg 15375 Ma^Mg^ 5565M2^M4 
256J2Mi^° ~ 256 ^ 512 ^ 256 

56O5M2M3M4 72875 Ms^ 3213 Ms 2515 M3M5 
256 J2 M/ 21504 J2 m/ ~ 256 J2 M/ ^ 512 J2 

21245 5929M2M6 5OO5M7 

^9216 J2 Mi^ ^ 1024 J2 ~ 3072 J2 ^' ^ 



^(ATS) _ 21023793 M2^ 12829887 M2^ M3 98342775 M2^ Mg^ 
^ ~ ~10240J3Mi^^ ^ 1024J3Mi^^ iogCkfW^^ 

4456305 M2'^M4 16200375 M2^ Mg^ 26413065 M2^ M3 M4 
512J3Mi^3 + i024J3Mi^2 + i024J3Mi^2 

12093543 M2^ Ms 83895625 M2 Mg^ 68294625 M2^ Mg^ M4 



2048J3Mi^2 32768 J3Mi^^ 4096 JSMi^^ 

12367845 M2^M4^ 13024935 M2^ M3 M5 15411627 M2^ Me 
2048J3Mi^^ 1024J3Mi^i 4096J3Mi^^ 



29 



32418925 Mg^ M4 17562825 M2 M3 M^^ 578655 M2 M^^ M5 

24:57 Q d^^Mi^^ ^ mOd^M^ ^ USd^M^ 
10050831 M2^M4 Ms 5472621 M3 Me 44207163 Ma^ M^ 
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2048J3M/° 1024J3M/° 20480 Mi 

151 1055 MaMs^ 7503125 M4^ 2642325 M3 M4 M5 



2048 J3 M/ 36864 M/ 2048 M/ 

11532675 Ms^ Me 6242775 M2 M4 Me 6968247 M2 M3 M^ 

16384 J3M/ 4096 J3M/ 4096 JSM/ 

4297293 M2^ Mg 12677665 M2 Mg 8437275 M5 Me 

4096 J3M/ ^ 32768 J3M/ ^ 32768 JSM/ 

8913905 M4M7 10156575 M3 Ms 8083075 Mio 

32768 J3M/ ^ 32768 J3M/ 98304 J3M/ ' 



(3.100) 



We note that the terms hsted above are only potentially relevant. For a m'th multi- 
critical model all terms involving M^, k > m will be subdominant in the double scaling 
limit and can be ignored. Obviously the arguments and results presented above can be 
translated to the case where the critical behaviour of the matrix model is associated 
with the endpoint, j/, of the eigenvalue distribution just be replacing M by J and dc 
by -dc. 

Let us turn now to the case of a symmetrical potential. In view of the scaling 
relations (3.95) it seems unnatural to work with terms like [p — x) and [p — y) and 
as shown in reference [47] there exists a way of avoiding this. However, an interesting 
relation can be derived if one analyses the scaling behaviour of Wg{p) in the x, y 
formalism. The details of this scaling analysis can be found in reference [47]. Here 
we just state the result. It turns out that in the case of a symmetrical potential the 
loop equation decouples completely into two independent equations. Each of these 
is a double scaled version of the loop equation for the asymmetric potential. One 
corresponds to the case where the critical behaviour is associated with the endpoint, 
x, of the eigenvalue distribution, the other to the case where the critical behaviour 
is associated with the endpoint, y. In particular the free energy of the symmetrical 
model can be written as a sum of two contributions one from each endpoint of the 
support of the eigenvalue distribution. The contribution from the endpoint, x takes the 
form of equation (3.75) and the contribution from the endpoint y appears from (3.75) 
when M is replaced by J and dc by —dc- However, since for a symmetrical potential 



A similar conclusion was reached in another approach in reference [58]. 

3.3 From 1-Matrix Model to (almost) r-function 
3.3.1 To Square or not to Square 

It is a well established fact that the partition function of the hermitian 1-matrix model 
in the double scaling limit is closely related to a r-function of the kdV hierarchy. In 



Jfc = {-i f+^Mk we find 






30 



references [59, 60] a continuum version of the loop equations of the symmetric hermitian 
matrix model was derived and it was found that the corresponding continuum partition 
function was the square of a r-function of the kdV hierarchy. A similar conclusion was 
reached in reference [61] where discrete loop equations were used as the starting point — 
but it appeared that just taking the double scaling limit of the partition function of 
the symmetric hermitian matrix model does not leave one with the square of a r- 
function. In an attempt to overcome these difficulties a renomalization procedure was 
prescribed. Equation (3.101) indicates that the reason why it is the r-function squared 
which appears in the here mentioned analyses is that a symmetrical potential for the 
matrix model is assumed. For a generic potential the continuum partition function 
should be related to the r-function itself rather than the r-function squared. As we 
shall see this can indeed be proven to be the case. 

After the analyses of references [59], [60] and [61] new development has taken 
place. A matrix model realization of a r-function of the kdV hierarchy — the Kont- 
sevich model — has been found [62]. Since this discovery there have been attempts 
to relate the partition function of the Kontsevich model to the partition function of 
various 1-matrix models. In reference [63] a limiting procedure which allows one, on 
path integral level, to pass from the partition function of the reduced hermitian matrix 
model to the square of the partition function of the Kontsevich model was presented. 
In reference [47] a different limiting procedure provided a way of passing from the par- 
tition function of the generic hermitian 1-matrix model to the non-squared partition 
function of the Kontsevich model. However, the latter limiting procedure was rather 
unconventional involving an analytic continuation of the size of the matrices entering 
the original 1-matrix model from N to —^N. Here we will show using the conven- 
tional double scaling prescription that it is possible to define continuum time variables 
for the generic hermitian 1-matrix model such that its correlators after appropriate 
renormalization — in compliance with the analysis of section 3.2.2 — turn into the cor- 
relators of the Kontsevich model expressed in terms of kdV times. Furthermore it 
appears that the double scaling limit of the partition function of the generic hermitian 
1-matrix model agrees with the partition function of the Kontsevich model (and not 
the square of the partition function) except for some complications at genus zero. The 
complications encountered for genus zero are of the same type as those encountered in 
reference [61]. However, we see no sensible way of removing these complications. The 
analysis which will be carried out below is based on reference [64]. 

3.3.2 The Kontsevich Model 

The Kontsevich model is defined by the partition function 
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where the integration is over N X N hermitian matrices. This partition function only 
depends on the parameters tk 



tfc = ^ Tr M-(2'=+^) (3.103) 

and expressed in terms of these it is a r-function of the kdV hierarchy [62] . As is the case 
for 1-matrix models the free energy has a genus expansion (cf. to equation (3.4)). 

The genus zero contribution to reads 

^6 2 iV2 ^.^^ \ m, + m, J ^ ' ^ 

where the m^'s are the eigenvalues of the matrix M and the parameter Uq is given by 
the boundary condition 

uo + ^E^=^==0. (3.105) 

It can be derived by means of the Dyson Schwinger equations of the model as done 
in references [65, 66, 67]. Alternatively it can be found by exploiting the fact that 
^Aorat|-jY^ {tfc}] is a T-function of the kdV hierarchy [68]. The higher genera contribu- 
tions can be written in the following form [67, 68] 

pKont ^ __L in(/^ - 1) (3.106) 

i^f = E («i • • -".klf " rrirh^ a > 2 (3.107) 

where the sum is over sets of indices obeying the following restrictions 

s 

J2[aj - 1) = 3g - 3, [a - s) = 2g - 2 (3.108) 

and where the moments Ik are defined by 

1 ^ 1 

^'=ivg K-2.„)>../. - 

The quantities (cti . . .as\a)^^°"'^ have a geometrical interpretation in terms of inter- 
section indices on the moduli space Aig^s of Riemann surfaces with g handles and s 
marked points [62, 69]. Different types of singular behaviour are possible for (3.107). 
These different types of singular behaviour can be indexed by an integer, m. As shown 
in reference [68] to reach the m'th type of critical behaviour one should neglect Ik with 
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k > keep constant, send /i — 1, /2, . . . , /m-i to zero and introduce a scaling 
parameter, s, such that 



remains finite while 



jq- 



l-h 



3/5 3/5 



3<q<m (3.110) 
■s (3.111) 



tends to zero. In this limit one has 



Fg = 1^ («1 -^^y </>2. (3.112) 

We see that this prescription for fine tuning the /'s is completely equivalent to the 
fine tuning of the Af 's in the vicinity of a m'th multi-critical point, playing the role 
of a(2™+i) (cf. to equation (3.86) and (3.87)). H owever in the case of the hermitian 
1-matrix model our starting point is the expression (3.38) for the free energy which as 
opposed to (3.107) contains terms which will never be of importance for any type of 
critical behaviour. Our conjecture is that with an appropriate definition of continuum 
times for the hermitian 1-matrix model the potentially relevant terms of (3.38) are 
exactly given by (3.107). 

3.3.3 Continuum Time Variables, Correlators and Loop Equations 

To motivate our choice of continuum time variables Tk for the non symmetrical her- 
mitian 1-matrix model let us write the boundary equation for the Kontsevich model 
as 

oo 

Y.ck{h + 8k,i){2uaf = Q (3.113) 

fc=0 

where was defined in (3.80). 

The idea is now to define Tk in such a way that by taking the double scaling limit 
of the boundary equations (3.19) and (3.20) we reproduce equation (3.113) with the 
Tfc's replacing the t^'s. Here and in the following we will consider a 1-matrix model 
for which the eigenvalues, at the critical point, are confined to the interval [j/c,a;c]. As 
mentioned earlier, when we move away from a given m'th multi-critical point by a 
change of coupling constants 8gi ~ (9(a™) in general both x and y will change. To 
keep the presentation as simple as possible we restrict ourselves to considering only the 
subclass of deformations for which y is kept fixed at j/c Expanding equation (3.19) 
in powers of [x — Xc) we find 

oo 

Y,c,{x-x,YM;{{g,])={) (3.114) 

p=0 



■*We note that this subclass of deformations does not include the deformations that lead to the 
emergence of the matrix model string equations (cf. to page 25). However the arguments of this and 
all following sections can be generalized to the case where y is not kept fixed. The expressions just 
become more involved. 
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where 

M;{{g,})=M,{x,,y,,{g,}). (3.115) 

Note that the couphng constants entering M!p{{gi]) are here completely arbitrary. 
Rewriting the boundary equation (3.20) we find 

oo 

Y: c,{x - x^y {x, M; {{g,}) + M;_, {{g,})} = 2. (3.116) 

Since Mp[{gi}) ~ a™"'', for < k < m the term Mp_^[{gi}) is subleading when 
compared to Mp[{gi}) except for M^^[{gi}) which is equal to 2 up to subleading terms 
of 0(a"+i). If we set 

X — Xc = a (2mo) (3.117) 
and define our continuum time variables by 

Tk + 4,1 = d'J' Ml {{g,}) , A; > (3.118) 

both (3.114) and (3.116) turn into the boundary equation of the Kontsevich model (3.113) 
(The reason why we include an additional factor ^/a d^J^ will become clear later.) 
We define continuum correlators W'^°'^^{tt\, . . . , tt^) by 

1 1 / 1 oo \ 

W{p)--V'{p) '-^ -hl/-*(,r)--^(rfc + 4,i)7r'=-'/' (3.119) 

2 « V ^k=o J 

W{pr,P.) + l^^—^ ^ ^f^-*(-i,vr,) + -i^^^,^) (3.120) 

W{p^,...,p,) '-^ ll1/-*(7ri,...,7r,) (3.121) 

a* 

where the scaling prescription (3.88) for pi is understood. Multiplicative renormaliza- 
tion is prescribed for all correlators. Additional subtractions appear in the case of the 
1- and 2-loop correlators. This is exactly the type of renormalization that we found to 
be necessary in the analysis of section 3.2.2 We note that the subtractions appearing 
above only concern the genus zero contributions to the 1- and 2-loop correlators. For 
the 2-loop correlator in the case where the two momenta coincide we use the following 
renormalization prescription 

Wip,p)'-^^^{w-\.,.) + ^). (3.122) 

The rationale for the multiplicative renormalization can also be given without referring 
explicitly to the detailed analysis of section 3.2.2. By expanding the moment Mp in 
powers of X — Xc one gets 



oo 



Mp = j: a-'-'" {Ti + 4,0 (2^o)'-^ r(p + i/2) (/ -p)! ' ^ " ^^'^^^^ 
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Using equation (3.123) the scaling relation (3.117) and (3.88) for x and p respectively 
and the boundary equation expressed in terms of T^'s one can by means of the chain 
rule show that the following relation holds 
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d.s.l. 



1 



d 



1 



dV{p) 



dV 



conti 



TT 



E(^ + l/2) 

fc=0 



1 
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7r'=+3/2 dn 



(3.124) 



Bearing in mind the relation (3.22) it appears natural to extract one power of for 
each loop in a given correlator. The relation (3.124) will be essential for the proof 
of our conjecture. (We note that from the equation (3.123) one can read off the 
continuum scaling behaviour of the moments for a given m'th multi-critical model 
and the relation (3.69) is easily reproduced.) Due to the peculiarities of the genus 
zero contributions to the 1- and 2-loop correlators it is convenient to use the loop 
equations in the genus expanded version (3.23). Let us introduce in equation (3.23) 
the continuum correlators. First we note that since Wo{p) itself does not scale we have 
in the double scaling limit 



K-2WMW,{p) 
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duj V'{uj 
c 27rz p — id 
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cont I 
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cont I 



TT — id 



(3.125) 



w 



here denotes a contour integral where the contour encircles infinity and where we 



have used the definition (3.119). To proceed let us write VK^f°"*(c<;) as 



W. 



cont I 



oo 



CO] = 2_^i_^-'l-3/2'^YCont,q_ 

q = 



(3.126) 



That W^°^^ allows an expansion of this type is obvious for g > I since for g > I we 



lave 



cont I 



dv^°^t(uj)^9^ where ^ycont(^) 



is given by equation (3.124). (It is also 



evident from the explicit expression for Wg{p) given in equation (3.51).) That the same 
is true for = will become clear in section (3.3.5). Performing the contour integral 
in (3.125) and making use of the definition (3.120) one obtains the following continuum 
loop equation 



2Wo=-*(vr)-^(rfc + 4,i)7r'="'/' 



^ 11/-*(7r)11/-;(7r) + 11/;!r(7r,7r) + 8,,, 

g' = l 



g,a 

1 

16^' 



(3.127) 
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3.3.4 Loop Equations and Correlators for the Kontsevich Model 

Inspired by the equation (3.124) let us introduce a loop insertion operator for the 
Kontsevich model by 



d 



1 



d 



and multi-loop correlators by 
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TTi, . . . , TTs 
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s> 1. 



7^1, 



(3.128) 



(3.129) 



Using the relation (3.103) it is easy to show by means of the chain rule that 
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(iyKont(^^ 

In this form the loop insertion operator can readily be applied to (3.104) to yield 



(3.130) 
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Wo'''°'^*(7ri,7r2) 



1 j TTi + 7r2 - 4mo TTi + 7r2 

4(7ri - 7r2)2 \ yvTi - 2MoV7r2 - 2uo ^rT^ 
Furthermore taking the limit tti — )• 7r2 in (3.132) we find 
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(3.131) 



(3.132) 



(3.133) 



4(7r - 2Mo)vr2 

We note that the 1-loop correlator of the Kontsevich model is analytic in the complex 
plane except for a square root branch cut [— oo,2mo] on the real axis. This is exactly 
the same analyticity structure as one obtains for the 1-loop correlator of the hermitian 
1-matrix model by the scaling prescriptions (3.117) and (3.88). With the definitions 
given above the master equation of the Kontsevich model can be written as 
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where 
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(3.134) 



(3.135) 



To make contact with the previous section we will rewrite this equation in a genus 
expanded version. To deal with the sum appearing on the left hand side of equa- 
tion (3.134) let us note that the genus g contribution to the 1-loop correlator can be 
expanded in the following way 
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Kont 
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(3.136) 
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Then making use of the definition of the t^'s, (3.103), we arrive at the foUowing form 
of the loop equation. 



w. 



Kont / 
g ' 



a-1 



g' = l 



(3.137) 



q,a 
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3.3.5 Proof of our Conjecture 

The task of proving statement made in section (3.3.1) amounts to proving the foUowing 
two identities 

Wo=°-*(7r) = Wo^"'°"*(7r), (3.138) 
Wr\T^i,T^2) = Wt'^\7r^,7r2) (3.139) 

where it is understood that the quantities on the left hand side should be expressed 
in terms of T^'s whereas those on the right hand side should be expressed in terms of 
tfc's. By comparing equation (3.127) and (3.137) it is easily seen that once this task 
has been fulfilled it follows by induction that 

W^^-'iTT, {Tk}) = W^'^-'iTT, {h}), g > 1 (3.140) 

since now the two sets of loop equations and corresponding boundary equations only 
differ by {Tk} appearing in one case and {tk} in the other. Furthermore by taking a 
glance at equation (3.124) and (3.128) bearing in mind the relations (3.22) and (3.119) 
- (3.121) one easily convinces oneself that 

11/-*(7ri, . . . , TT,, {n}) = Wf°-\T:r, . . . , TT,, {h}), g,s>l. (3.141) 

From equation (3.119) and (3.140) it follows that 

W,{p) '-^ -Wf°^\T^), g > 1 (3.142) 
a 

since both the second term on the left hand side and the second term on the right hand 
side of (3.119) are of zeroth order in genus. Now due to the similarity between ^yeott(^) 
and ^yicf„t(^) we immediately find 

To address the g = case we note that equation (3.120) and (3.139) imply that the 
double scaling limit of Wo{pi,P2) differs from W|j^°"*(7ri, 'K2) only by a term which does 
not depend on any couplings. Therefore we have 

Wr\T^i, • • • , vr.) = Wf^'iTT^, . . . , TT,), 5 > 3. (3.144) 
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However, equation (3.138) does not allow us to conclude anything about the relation 
between the genus zero contributions to the partition functions of the two models 
because of the subtractions appearing in equation (3.119). 

Let us now turn to the proof of the relations (3.138) and (3.139). The proof of the 
latter is by far the most straightforward since Wo(pi,P2) is universal, i.e. it does not 
contain any explicit reference to the coupling constants. Taking the double scaling limit 
of (3.55) following the prescriptions (3.117) and (3.88) one easily reproduces the genus 
zero part of (3.120) with Vl^o^°"*(7ri, 7r2) replacing Wo^°"*(7ri, 7r2). Likewise taking the 
double scaling limit of (3.56) one gets exactly the genus zero part of the relation (3.122) 
with Wq°"'^ {tt , tt) replaced by W|j^°"*(7r, tt). To prove the relation (3.138) we will prove 
that 

1 11 oo 1 

Wo{p)-^V'{p) + --J2{T, + 6,,,)7r'-'/' '-^ -Wt'^'iTT). (3.145) 

2 a 2 a 

First we note that due to equation (3.18) we can write the two first terms as 

Wo{p) - Iv'ip) = f { I'u". 'A \ ■ (3-146) 

Furthermore inserting our definition of continuum times (3.118) into the remaining 
term on the left hand side of (3.145) we get 
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^ 2 fc=o 2 fc=o 

Wo{p) - Iv'ip)] + \{p - x^Y'l^ip - y^y/'M^ {{g,}) 
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x=xc,y=yc 

1/2 



2 Joo 27rz p — (jj [{p — Xc){io — yc) J 

To obtain the second equality sign we have made use of the rewriting of Wo{p) given 
in (3.27) and the scaling relation for p, (3.88). (In particular we have used that [p — y) 
in the continuum limit can be replaced by dc ) We note that it was in order to be 
able to carry out this step that we had to multiply our boundary equation by d^J^^fa. 
To obtain the third equality sign we have made use of the relation (3.146). So our 
statement is now the following 



Wo^^'"-^* = lim a • 

d.s.l 



1 ( dio V'{uj) i \I{P- x){p - yc) ^(cu - Xc){p - yc) 



(3.148) 

The similarity with equation (3.131) is striking and the equality is straightforward to 
prove. To do so one expands {p—uj)~^ in powers of (^5^) ^^'^ quantities (p—x) and 



[(jj — x) in powers of (^^zf^j and y ^J^^ j respectively. The factor [p — yc) can simply 
be replaced by dc- The ^ term of (3.148) vanishes as it should. This is actually 
ensured by the boundary equation. In the process of expanding the integrand it proves 
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convenient to pull out a factor {p — x)^/^. The result of the expansion procedure is the 
following expression for the right hand side of (3.148) 



(. 6=1 m=l v" ) 

= li--^^of'f:'t^^^^AT^ + ^rn^). (3.150) 
6=1 m=l 

By rewriting (3.131) using the definition (3.103) of the time variables it is easy to show 
that the expression (3.150) is exactly Wq'^°"'^{'k). 

Now comparing the expression for Afp, (3.123), with the expression for /p, (3.109), 
bearing in the mind the relation (3.103) we see that (cti . . . = (cti . . . as\a)g 

and hence that we can carry over the geometrical interpretation of the former to the 
latter. 



3.3.6 Comparison with other Approaches 

Although it is already clear from equation (3.101) let us take the opportunity to show 
explicitly that there is no contradiction between our result and the results of refer- 
ences [59, 60, 61, 63]. The strategy applied to the generic hermitian matrix model in 
sections 3.3.3-3.3.5 allows us to address the symmetrical model as well. However one 
has to use the complex matrix model as an intermediate step. This is of course due to 
the fact that the symmetrical hermitian matrix model does not contain a complete set 
of operators. The complex matrix model does — and in reference [47] it was shown that 
in the double scaling limit the partition function of the complex matrix model involv- 
ing matrices of size N X N equals the partition function of the symmetrical hermitian 
matrix model involving matrices of size 2N X 2N. Namely, it holds that 

= 4^^f ^ (^--M (3-151) 

where C refers to the complex matrix model. The complex matrix model is in many 
respects very similar to the generic hermitian matrix model. It has a set of loop equa- 
tions which can be written in the same form as that of equation (3.15). The appropriate 
requirement concerning the analyticity structure of its 1-loop correlator is that it has 
only one square root branch cut \—\fz^\fz\ on the real axis. This corresponds to the 
eigenvalues of the matrix (^V) being confined to the interval [0, z\. With this require- 
ment one can solve the loop equations genus by genus. The solution of course depends 
on the parameter z which is determined by a boundary condition similar to (3.20). As 
before the multi-loop correlators can be found by applying a loop insertion operator 
to the 1-loop correlator and as before expressing the higher genera contributions to 
the correlators is facilitated by introducing a moment description. The m'th multi- 
critical point is reached when the distribution of eigenvalues of (^V) acquires (m — 1) 
extra zeros at the endpoint z. To relate the double scaling of the partition function 
of the complex matrix model to the one of the Kontsevich model one takes the same 
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line of action as for the generic hermitian matrix model. Appropriate continuum time 
variables are defined by the requirement that the boundary equation of the complex 
matrix model reproduces the boundary equation of the Kontsevich model when the 
double scaling limit is taken. The resulting time variables turn out to be related to the 
moments of the complex matrix model by an equation similar to (3.f f 8). For the corre- 
lators the appropriate renormalization prescriptions are very similar to (3.f f 9)-(3.f 2f ) 
and furthermore for the loop insertion operator one has again a relation like (3.f24). 
However, a closer analysis of the loop equations shows that in stead of (3.f 43) we have 



By comparing (3.f52) with (3.f5f) we immediately see that except for some compli- 
cations at genus zero we get in the double scaling limit from the partition function of 
the symmetrical hermitian matrix model the square of the partition function of the 
Kontsevich model. A similar statement of course holds for the complex matrix model, 
the only difference being that starting from complex matrices of size N X N one ends 
up with the Kontsevich model involving matrices of size 2N X 2N. 

Let us close this section by briefly discussing the unconventional limiting procedure 
mentioned in section (3.3.1). This procedure consists in analytically continuing the size 
of the matrices of the generic hermitian 1-matrix model from N to ■ N = —j^N, 
rescaling matrices and couplings appropriately and sending e to zero. With this pre- 
scription one can turn the path integral defining the generic hermitian 1-matrix model 
into the path integral defining the Kontsevich model. However, solving the boundary 
equations (3.19) and (3.20) in the limit e — )• assuming a non symmetric support of 
the eigenvalue distribution also allows us to apply the e — )• procedure to the represen- 
tation (3.38) of the free energy of the generic hermitian 1-matrix model. Following this 
line of action one finds that only terms of (3.38) which fulfill the requirements (3.72) 
survive in the limit e — )• 0. Hence the e — )• procedure leaves us with exactly the 
terms which are potentially relevant for the double scaling limit! In addition the lim- 
iting procedure sends a given of the surviving terms of (3.38) into the term of (3.107) 
which appears from it when (Af^ -|- Sk^i) is replaced by Ik and dc removed. We refer 
to [47] for details. The e — )• prescription is an amazing limiting procedure. It picks 
out exactly the terms which are relevant for the double scaling limit and yet we can 
not relate e in any way to the double scaling parameter, a. In addition it allows us 
to access less explored regions of the moduli spaces Aig^s of Riemann surfaces with g 
handles and s marked points. In the same way as the model arising in the limit e — )• 0, 
i.e. the Kontsevich model, can be viewed as a generating functional for intersection 
indices on moduli spaces A^p,s, the model corresponding to a finite e can be viewed as 
a generating functional for intersection indices on discretized versions of moduli spaces, 
e being the step of discretization. This gives rise to a geometrical interpretation of the 
coefficients of those terms in (3.38) which are subdominant in the limit e — )• and 
hence in the double scaling limit. For a discretized moduli space one must take into 
consideration contributions from the boundary when calculating intersection indices. 
The boundary of the moduli space Aig^s appears by the so-called Deligne-Mumford 
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compactification [70] and consists of (products of) lower dimensional moduli spaces. 
While the coefficients of the leading order terms of Fg are related to intersection indices 
on moduli spaces A^p,s, s < 3(/ — 3, the coefficients of the subleading terms are related 
to intersection indices on the boundary of these spaces. For a description of the idea 
of the discretized moduli space we refer to [71, 72]. 

3.4 A Remark on the Two-Matrix Model 

The two-matrix model has turned out to encode a vast amount of information about 
the interaction of matter fields with 2-dimensional quantum gravity. This property is 
most clearly exposed by the following two-matrix integral [28] 

Z[6,A,iV] = e^[^'^'^] 

= f d<P+d<P_exp\^-NTT (^^^l + ^^^_-b<P+<P_-X[<Pl + <P^_)'^Y (3.153) 

By the perturbative expansion of F[b^X^N] one generates the same simplicial mani- 
folds as by the perturbative expansion of F[A,A^] appearing in (3.1). However each 
triangle is now equipped with a label which is either -|- or — , and the weight of a given 
triangulation is modified by a factor ^1:^2 foi" each pair of neighbouring triangles with 
identical labels and a factor for each pair of neighbouring triangles with different 
labels. This clearly corresponds to a situation where we have added to the Einstein 
Hilbert action (2.12) the Ising model action (2.19) with 

e-'f = b, (3.154) 

the spins being placed in the centres of the triangles. The genus zero contribution to 
the free energy of the model (3.153), Fo[^, ^] can be calculated using the technique 
of orthogonal polynomials [73, 74]. Analyzing the singularity structure of Fo[b^ A] one 
finds that there exists a point in the coupling constant plane for which the spin system 
has a third order phase transition and the geometrical system becomes critical in the 
sense described in section 2.3 [28, 29]. In the vicinity of this point the model describes 
unitary conformal matter with c = | coupled to 2-dimensional quantum gravity. 
The generic hermitian two-matrix model takes the following form 

Z[c,{g,},{g,},N] = e^W^'>'«'>'^] 

= I dXdY exp {-N Tr [V(X) + V(Y) - cXY) } (3.155) 

where 

V(X) =J2^X^, V(Y) = J2^Y'. (3.156) 

j=i J j=i J 

Recently it has been shown that in the coupling constant space of this model one can 
find critical points describing the interaction of any rational conformal matter field 
with quantum gravity [34]. For rational (p^q) matter the critical potentials Vc and Vc 
can be taken to be of degree p and q respectively. For unitary conformal matter with 
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(p, q) = (m, m + 1) one also has the possibihty of choosing a critical point for which the 
corresponding potentials are identical and of order m. Hence all the unitary models 
can be viewed as Ising models living on random surfaces. 

In the light of this recent development it would be interesting to try to generalize 
the techniques described in the previous sections to the two-matrix case. As mentioned 
above the genus zero contribution to the free energy for the model given by (3.153) 
has been calculated using the method of orthogonal polynomials. The study of loop 
equations for this model has also been initiated [75, 76, 77] but explicit results for 
correlators are few and limited to genus zero. We are of the opinion that a more 
detailed analysis along the lines of the previous sections should indeed be possible and 
might lead to a better understanding of the structure of two-matrix models both away 
from and in the continuum. Let us try to describe our present understanding of how 
such an analysis should be carried out [79]. The starting point would of course be 
the loop equations of the generic hermitian two-matrix model. To derive these in an 
appropriate form it is convenient to consider the following redefinition of the field, X 

X^X + e^-^. (3.157) 

p — X q — y 

Under a transformation of this type the measure changes as 



dX ^ dxll+eTi {^—-] Tt {^—-^—-]\ (3.158) 



and we get to the first order in e 

dXdYlTi (^^) Tr (^^^tt) -iVTr {V'{X)- ^ 




p — X q — Y J \ p — X q — Y J 

+cNTt (f^— Uexp{-iVTr {v{X) + V{Y) - cXY)} . (3.159) 



Gpg = {TT- ^^^) (3.160) 



For the two-matrix model we have of course a new type of correlators, namely correla- 
tors involving both X and Y fields. Let us introduce the following definitions 

1 1 

p^^Y' 

and 

G'™(p) = (Tr F"), (?"(?) = (TrX"^—). (3.161) 

p — X q — Y 

Then it obviously holds that 

G„ = f = f ,3.162) 

ri=0 " m = " 

Furthermore let us define loop insertion operators by 

dV{p) f^.P^+^dg,' dV{q) fr{q^+'dg/ > 
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Assuming that eigenvalue distributions for X and Y have been introduced we can now 
write the equation (3.159) as 

^ ^^G., = G\p)G,, + c{qG,, - G\p)) + -^,G,, (3.164) 

C ZTTZ p — (jj dV [p) 

where C is a curve which encircles the supports of the eigenvalue distributions and 
where we have made use of a factorization condition similar to (3.16). As in equa- 
tion (3.15) the term -j^fj^^Gpq is suppressed by a factor when compared to the rest. 
Let us describe how one can in principle solve equation (3.164) step by step. First we 
take the coefficient of ^ in the expansion of (3.164) in powers of ^. This gives 

oo 

Y,g,G,.r{q)=c{qG%q)-l). (3.165) 

Of course (3.164) has a twin equation where the roles of X and Y are interchanged. 
Expanding this twin equation in powers of ^ and taking the coefficient of ^ we find 

the following recursion relation for G"'{q) 

cG-{q) = I ^^G-^LO) - G'{q)G-\q) - -^G^^q). (3.166) 
Jc ZTTi q - u dV{q) 

By means of this relation we can express G"'{q) entirely in terms of G^{q). Neglect- 
ing the subleading term in (3.166) we that see that this equation allows us to obtain 
from (3.165) an algebraic equation for the genus zero contribution to the 1-loop corre- 
lator for the field Y. The degree of the algebraic equation is the same as the degree 
of the potential V{X). This is in full compliance with the situation in the 1-matrix 
case where the genus zero contribution to the 1-loop correlator was determined by an 
algebraic equation of the second degree since a 1-matrix model can be considered a 
two matrix model where one potential is gaussian and has been integrated out. The 
equation that determines G^{q) contains a number of unknown integration constants. 
These constants must be determined by imposing appropriate requirements on the 
analyticity structure of G^{q). A similar situation was encountered in the 1-matrix 
case. Of course, unless special circumstances occur we can only solve the equation for 
G^{q) when the degree of the potential V{X) is less than or equal to four. However as 
long as the degree of V{X) does not exceed four we can probably find and expression 
analogous to (3.18) giving G^{q) for arbitrary V{Y). In case we succeed in arriving at 
such a formula we have the possibility of relating the various types of critical behaviour 
of the 2-matrix model to the behaviour of its eigenvalue distributions in analogy with 
what was done for the 1-matrix model. Furthermore such a formula might provide us 
with an important clue as to what are the fundamental variables of the two-matrix 
model. Maybe a set of generalized moments will reveal themselves when the 1-loop 
correlator is analysed in detail. This was what happened in the 1-matrix case (cf. to 
equation (3.27)). Finally a formula for G^{q) for generic V{Y) might also allow us 
to develop an iterative procedure for calculating higher genera contributions to G^{q) 
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similar to the one described for the 1-matrix modeL We can derive a genus expanded 
version of the equation for G^{q) combining (3.165) and (3.166). With a generic V{Y) 
we have a loop insertion operator ^^^^^ and the equation closes. Probably the iterative 
solution is again conveniently expressed in terms of the eigenvectors of some linear 
operator. When we have calculated the genus zero contribution to G^{q) we can in 
principle determine Gpq from the twin equation of equation (3.164). While this is 
certainly doable when the potential V{Y) has a finite number of terms it might be 
practically impossible for a generic V{Y). However, if it is not, an iterative procedure 
for calculating higher genera contributions to Gpq is very likely to exist since the struc- 
ture of (3.164) is very similar to that of (3.15). From G^{q) and Gpq a whole series of 
correlators can be calculated by application of the loop insertion operator ^J^^ ^ and the 

operator ^. The s-loop correlators for the X field seem difficult to access, however. 

To our opinion one should not be deterred by the seemingly large degree of com- 
plexity of the loop equations of the two-matrix model. It might very well happen that 
a relatively simple structure of correlators and the free energy emerges as it was seen 
in the case of the 1-matrix model. Studying the loop equations of the two matrix 
model might lead to new insight concerning rational conformal theories coupled to 2- 
dimensional quantum gravity. Furthermore the same loop equations have been found 
to turn up in the Kazakov-Migdal model of induced QCD [78]. 



3.5 Non Perturbative 2D Quantum Gravity ? 
3.5.1 Non Perturbative Effects in Matrix Models 

The wish to understand in detail the the simple model (3.1) has brought us far. We have 
found the complete perturbative solution of the generic hermitian 1-matrix model (3.6), 
i.e. we know (at least in principle) all terms in the genus expansion of its free energy 
as well as of any of its correlators. We have localized and classified all critical points 
in the coupling constant space of the generic model and devised a prescription, namely 
the double scaling prescription for how to obtain continuum theories in the vicinity 
of these points. In the same sense as for the discrete models we have determined the 
complete perturbative solution of these continuum theories, i.e. we have access to (in 
principle) all terms in the genus or double scaling expansion of their free energy and 
their correlators. However, as regards the understanding of 2-dimensional quantum 
gravity beyond the topological expansion our achievements are less impressive. Most 
efforts have naturally been concentrated on studying the pure gravity, c = 0, case. Let 
us briefiy summarize the results of these efforts. For convenience we will refer to a 
specific m = 2 multi-critical point, namely the one given by 

9i = ^s^ ^3 = "3' ^/.' = 0, for z>3, z = 2. (3.167) 
This critical point can be reached starting from the following matrix model potential 

¥{<!>)= g<l>-^<l>' (3.168) 
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which appears from the one in (3.1) by a simple rescaling. Each term in the topological 
expansion of the model with the potential (3.168) is a power series in which 
is convergent when g > (/c, divergent otherwise. Let us (for g > gc) set g — gc = a^^R- 
Then the double scaling limit is defined as the limit a — )• 0, — > oo, = a^N^ fixed. 
The topological expansion in this limit becomes an expansion in the parameter GA^^^^. 
All attempts to perform the sum over topologies for this expansion have failed. The 
series is not even Borel summable. (We note that the divergencies appear only when 
we try to sum over topologies and are entirely due to the rapid increase of the entropy 
of manifolds with genus.) The Borel non summability of the perturbation series for 
the model given by the potential (3.168) can be understood as a consequence of the 
action not being bounded from below and hence the original matrix integral being ill 
defined. (Let us mention that it is actually possible to reach a m = 2 multi-critical 
point starting from a potential which is bounded from below i.e. a potential of the 
type Yli=i gi^' with g2j > 0. We will not enter into a detailed discussion of these more 
complicated m = 2 multi-critical models. Let us just note that having a matrix model 
coupling constant gi > corresponds to considering space time manifolds where z-gons 
appear with negative weight.) For some time it was the hope that the connection 
between the double scaling limit of the 1-matrix model and the kdV hierarchy would 
lead to a non perturbative definition of 2D quantum gravity. For the model given by 
the potential (3.168) this connection implies that in the double scaling limit its specific 
heat, C{Ar) = —f^ obeys the following differential equation — the Painleve equation 



-^C" + C^ = Kr. (3.169) 
6 

However this observation does not lead to very much progress as regards our search 
for a consistent non perturbative definition of 2D quantum gravity. There is an infi- 
nite number of solutions to the Painleve equation which have the correct asymptotic 

1/2 

behaviour, namely C{A.b) — )• A^ for A^^ — )• 00 (cf. to equation (3.5)). The difference, 
(5C, between two such solutions is invisible in perturbation theory. Its asymptotic 
behaviour can be found by linearizing (3.169). The result reads 

8C ~ A-^/^xp (^-^G-i/^Af ^ . (3.170) 

Bearing in mind that our expansion parameter in the double scaling limit is GA^ 
it is obvious that the topological expansion can not provide us with any information 
which would allow us to choose between the different solutions. Additional physical 
constraints are needed. The search for such constraints has so far been unsuccessful. 
The above mentioned ambiguity is not the only problem, however. All real solutions 
to the Painleve equation have an infinite number of poles on the negative real axis. 
The existence of these poles makes it impossible to define continuum s-loop correlators 
without violating the continuum loop equations [53]. 

The argument of the exponential in (3.170) can be recognized as (minus) the action 
of a real instanton present in the original matrix model [53, 80]. The presence of this 
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real instanton explains the non Borel summability ol the perturbation series. In the 
large- limit the eigenvalues of the matrix (f> are confined to an interval [y^x] on the 
real axis. The instanton configuration is one where the largest eigenvalue has been 
moved from the point x to some point x > x while the others remain at their original 
positions. This change of configuration can be described by the following change in 
the eigenvalue distribution 

A (uo(A)) = ^ [S{X - 5) - S{X - x)] . (3.171) 

Under a deformation of Uo[X) of this type the effective action, S*, of the matrix model, 
i.e. —Fq changes (to the leading order in N) by 



AS 



N"^ j dXA (mo(A)) X (^y(A) - 2 y" djiuoiii) In |A - 
N £ dX (v'{X) - 2 J dfiUo{fi)j^^ 



N / dX{-2Wo{X) +V'{X)) 

-2N r dXW^'^'iX). (3.172) 



To arrive at the last equality sign we have once again made use of reference [14]. We 
note that Wq*^"'(A) is exactly the scaling part of Wo(A). For the potential (3.168) we 
immediately find using the relation (3.27) 

Wr\X) = -i(A - xyf\X - y fl\X + \{x + y)) (3.173) 

where x and y are given by (3.19) and (3.20). The instanton configuration referred to 
above is characterized by = '^^^'^^ oc Wq'^'^\x) and corresponds to the situation 

i = -\{x + y). (3.174) 

It is easily verified that x = —^{x-\-y) > x, that 5; is a maximum for AS* and that 
AS* — > — oo for X — )• oo. Hence there is a barrier of height AS{x = —\{x + y)) which 
prevents the eigenvalues from escaping to infinity. The probability that the largest 
eigenvalue tunnels through this barrier is given by exp( — AS'(5: = —^{x + y)). In the 
large- limit this probability is exponentially suppressed and so is the contribution to 
the free energy from the instanton configuration. However, in the double scaling limit 
the situation is different. In this limit one finds using the boundary equations (3.19) 
and (3.20) 



X — Xc = —2xc y j (3.175) 

y-yc = 0{x-x,y (3.176) 
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(3.177) 



\ 9c J 

and inserting (3.175), (3.176) and (3.177) in (3.172) leads to 




(3.178) 



{^S)d.s.l. 



N 




\ 9c J J-2 
1/2 a5/4 / 12^3^ 



\ 5/4 ^ 

] I dz^/z + 2 (z - 1) 



G 



(3.179) 



We see that in the double scahng hmit the action is influenced by the presence of 
the instanton and that as advertised the amount by which the action is modifled is 
exactly equal to minus the argument of the exponential entering the non perturbative 
ambiguity of the Painleve equation. As noted earlier the presence of the real instanton 
is what obstructs the Borel summation of the perturbation series. 

Analyses of the properties of the perturbation series have also been carried out 
for matrix models describing gravity coupled to various matter flelds [81]. For the 
m'th multi-critical models corresponding to gravity coupled to conformal matter with 
(p, q) = (2, 2m — 1) the perturbation series has been found to be Borel summable if m is 
odd and not Borel summable if m is even. The non Borel summability when m is even 
can in all cases be explained by the presence of a real instanton. The corresponding 
instanton conflguration is of the same type as described for m = 2 above. In the case 
of the m odd models there are no real instantons. For a discussion of the general 
{p^q) case we refer to [81]. Let us just mention that for all the unitary models the 
perturbation series is not Borel summable. 

In the case of pure gravity there have been several attempts to overcome the above 
mentioned problems by modifying the matrix model construction, the idea being that 
the modifled theory should give a non perturbative deflnition of 2D quantum gravity 
while containing the same perturbative information as the original theory. Since the 
problems with the model given by the potential (3.168) can be traced back to the 
fact that the corresponding matrix integral is ill deflned one way to proceed is to 
try to modify the theory directly on path integral level. In references [53, 80, 82] 
it was pointed out that one can obtain a well deflned matrix integral with the same 
perturbative expansion as the original one with the potential (3.168) by choosing as the 
contour for the integration over eigenvalues in stead of the real axis a path which goes 
from — oo to e'f This procedure of course also ensures that all correlators are well 
deflned but they now have a non vanishing imaginary part. This holds in particular 
for the speciflc heat of the model which can be shown to correspond to a complex 
solution of the Painleve equation. This solution is unique and has no poles on the 
real axis. However, due to the appearance of imaginary parts in the partition function 



47 



and the correlators the contour deformation strategy does not lead to an immediately 
physically acceptable theory of quantum gravity. 

Another strategy for constructing a non perturbative extension of the original model 
has been suggested by Dalley et al. [83, 84]. It consists in replacing the Painleve 
equation by a differential equation of order four. This differential equation can be 
derived along the lines of the original derivation of the Painleve equation for the model 
given by (3.6) and (3.168) using as the starting point a complex matrix model with a 
potential, V{<f>^<f>), of the third degree in {4>^4>). For the eigenvalues A? of {4>^4>) there 
is a lower cutoff, namely A^^ > 0. The critical point of interest is characterized by the 
eigenvalues of {4>^4>) being confined to the interval [0, z] and the eigenvalue distribution 
acquiring two zeros at the endpoint 0, i.e. at the cutoff, in addition to its usual square 
root singularity at the endpoint z. This is a critical point of which no analogue exists 
for the hermitian matrix model since for hermitian matrices there is no lower cutoff 
on the eigenvalues. The above mentioned fourth order differential equation can also 
be derived without referring to a specific model but using the properties of the kdV 
fiow equations [84]. This differential equation has solutions which have the correct 
asymptotic expansion for A^^ — )• oo. Numerical investigations indicate that there is 
only one real solution with this property and that this solution does not have any poles 
on the real axis. 

In the following we will describe yet another proposal for defining 2D quantum grav- 
ity non perturbatively which we will denote stochastic stabilization. It is based on the 
prescription of Greensite and Halpern for stabilizing bottomless action theories [85]. 
This prescription results in the original ill defined J-dimensional quantum field the- 
ory being replaced by a well defined [d -\- 1) dimensional quantum field theory which 
contains the same perturbative information as the original one. In the case of the zero- 
dimensional matrix model given by (3.6) and (3.168) the corresponding I-dimensional 
theory corresponds to a simple quantum mechanical system. Stochastic stabilization 
in the context of matrix models has been studied by various groups [86]-[98]. For the 
stabilized model the partition function as well as all correlators are automatically real. 
Hence the non perturbative definition of 2D quantum gravity provided by stochastic 
stabilization differs from the one provided by the contour deformation method. A 
detailed comparison of correlators obtained by the two methods was performed in ref- 
erence [90]. In reference [91] stochastically stabilized 2D gravity was compared to the 
non perturbative definition of 2D gravity advocated by Dalley et. al. and it was found 
that these two approaches do not agree either. In both approaches the original geo- 
metrical picture is lost. With stochastic stabilization we loose in addition the kdV fiow 
structure. However, stochastic stabilization gives rise to a simple quantum mechanical 
interpretation of several matrix model characteristica. We will mention some examples 
in the next section. In particular we will reinterpret the non perturbative ambiguity 
of the Painleve equation. Stochastic stabilization in addition seems to present to us 
the possibility of a strong coupling expansion of 2D quantum gravity. This will be the 
subject of section 3.5.5. 
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3.5.2 Stochastic Stabilization 



Stochastic stabilization is a prescription which assigns to a J- dimensional bottomless 
action theory a well defined [d + l)-dimensional quantum field theory with the same 
perturbative information as the original ill defined theory. It relies on the observation 
that the expectation value of an operator Q in a J-dimensional field theory with an 
action, S[(p]^ bounded from below can be viewed as the expectation value of the same 
operator in the ground state of a ( J + l)-dimensional quantum field theory. 

<Q>=^fd<P e-'^^y^Q[<P] = fd<P ^f'M Q[<P]. (3.180) 
The [d + l)-dimensional theory is described by the Fokker-Planck Hamiltonian 



H 



FP 



d'^X 




1 S^S 



(3.181; 



which is a positive semi-definite operator 



Hfp = I d'^x R\x)R{x), R = t-^+ ' 



2n 



It is easily verified that 



-S[<f>]/2h 



Z 



(3.182) 



(3.183) 



is a normalizable eigenvector of Hpp corresponding to the eigenvalue £"0 = and hence 
qualifies as the ground state of the [d + l)-dimensional theory. However, if the action 
S[(p] is not bounded from below the wave functional given by (3.183) is not normalizable 
and does not qualify as the ground state of Hpp. The Fokker-Planck Hamiltonian is 
nevertheless still a positive semi-definite operator. Hence it also in this case has a well 
defined ground state, $o(^), now corresponding to a vacuum energy Eq > 0. This 
ground state, $o(^), allows us to define a stabilized version of the original ill defined 
theory, namely the theory with action Seff given by 



MM 



-S,ff[,f>]/2h 



Z, 



(3.184) 



eff 



For the stabilized theory expectation values are defined by 

1 



(Q) 



z 



eff 



{MQ\M- 



(3.185) 



It can be shown that the stabilized theory has the same classical limit, the same 
perturbative expansion in coupling constants and the same ^^-expansion as the original 
theory. In addition the action of the stabilized theory is of course bounded from below. 
Studying the stabilized theory amounts to studying the [d -\- l)-dimensional theory 
described by the Hamiltonian Hpp. 
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It is obvious that stochastic stabihzation offers to us the possibility of a non per- 
turbative definition of 2D quantum gravity. One can easily convince oneself that for 
the zero- dimensional matrix model defined by (3.6) Hpp takes the following form 

H,, = NiR'R). «.. = >(i^ + ^£). (3.186, 

Since our vacuum must be U (N) invariant and since we will only be interested in U (N) 
invariant observables we can write (3.185) as 



N 

(g) = y n dh lii^^ - ^if $^({Afc})g({A,}) (3.187) 

k=l i<j 



where {Xi} are the eigenvalues of the matrix (p. Likewise we can diagonalize the operator 
Hpp. In this picture the wavefunction we are interested in, namely 

$o({Afc}) = n(A. - A,) $o({Afc}) (3.188) 

i<j 

is totally antisymmetric. Hence the 1-dimensional well defined theory which by the 
stochastic stabilization scheme gets associated with our original zero- dimensional ma- 
trix model in fact describes a Fermi gas consisting of N particles. Let us now specialize 
to the matrix model with the potential (3.168). In this case the Hamiltonian of the 
fermionic system reads 

N 

Hfp = Nj2Hf,[X,l (3.189) 

8 = 1 

1 

HfAM = -]^^ + V/.(A), (3.190) 

V/.(A) = lig-\'Y + \. (3.191) 

The situation is simple. The fermions do not interact. The ground state of the system 
is hence the Slater determinant of the N lowest eigenstates of the single particle Hamil- 
tonian Hfp and the vacuum energy the sum of the corresponding eigenvalues times N. 
For potentials of degree larger than three the operator J^^^ gives rise to interactions 
between eigenvalues and the situation is more complicated. 

The model given by (3.189) - (3.191) allows us to attribute a meaning to the par- 
tition function and the correlators of pure 2D gravity for any value of the coupling 
constant g and without making any restrictions on topology. Hence we have a truly 
non perturbative definition of 2D quantum gravity. We note that the stochastic stabi- 
lization scheme always results in real answers for correlators and the partition function. 
No exact solution is known for the quantum mechanical problem given by Hpp. In 
reference [89] a recipe for numerical calculations of correlators of 2D quantum gravity 
based on a numerical solution of the quantum mechanical problem was given and some 
explicit results obtained. However, there is also quite a lot which can be learnt by 
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studying the model given by (3.189)- (3.191) by semi-classical techniques. In particu- 
lar many ol the characteristica ol the underlying matrix model can be reproduced and 
given a simple quantum mechanical interpretation. We will present some examples in 
the following sections. 

3.5.3 The Large-iV limit 

Let us consider the large- limit of the fermionic system given by (3.189)-(3.191). 
From (3.190) it appears that plays the role of . Hence we see that the genus 
expansion of our original matrix model is nothing but the WKB expansion of the 
quantum mechanical system defined by (3.189)-(3.191). In reference [98] the effect of 
applying the WKB approximation to the quantum mechanical system was analysed in 
some detail. Here we will just mention the main results. Since we are dealing with a 
system consisting of N fermions let us denote the energy of the A^'th level of Hjp by 
Ep. Furthermore let us assume that for E < Ep there is only one set of turning points, 
{y{E), x[E))^ for a classical particle with energy E moving in the potential Vjp. In the 
WKB approximation (with h = Ep is determined by 

1 = - / dXJEp - V/p(A). (3.192) 

Now going back to the expression (3.29) one finds that for the matrix model with the 
potential (3.168) the spherical level eigenvalue distribution is given by 

t^o(A) = -Jc{x,y) - V/p(A), y<\<x, g > gc (3.193) 

TT ^ 

where c(a;, y) is equal to ( Tr (f)) and can by expressed in terms of x and y as 

y) = Tr ^) = Q(x - y)^ + i Q(x^ - y')^. (3.194) 

Bearing in mind that Uo[x) = Uo[y) = and Uo[X) > for A it follows from the 

normalization condition (3.17) that for g > gc, Ep can be identified with the constant 
c[x,y). This implies that 

Vjp{X) - Ep = [wr'Wy , g>gc, X>x (3.195) 

where we have made use of the relation between the eigenvalue distribution and the 
1-loop-correlator (cf. to equation (3.27) and (3.29)). From the expression (3.173) for 
Wq*^"'(A) we see that for g > gc we have the situation depicted in figure 2a. As already 
mentioned several times the critical point of the matrix model is characterized by the 
eigenvalue distribution Mo(A) acquiring an extra zero at the endpoint x. This means 
that as we let g approach gc from above the point —^{x-\-y) will approach the point 
X until at g = gc they merge into a single point Xc (cf. to (3.173)). We see that in 
the quantum mechanical picture the merging of the two points corresponds to the local 
minimum of the potential disappearing and the potential acquiring a horizontal tangent 
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Figure 2. The figure to tlie left (fig. 2a) represents tlie situation g > gc- 
In tlie other figure (fig. 2b) g < gc- The dashed lines represent the Fermi 
energy in the two cases. 

For g < gc we have the situation depicted in figure 2b. The potential now only has one 
stationary point, namely the point corresponding to the global minimum. It is evident 
that, as pointed out in the previous section, the quantum mechanical system still has 
a well defined ground state. As opposed to this our original matrix model makes 
no sense for g < gc- Fven the planar approximation breaks down. In the quantum 
mechanical picture this breakdown is signalled by the Fokker-Planck potential acquiring 
a horizontal tangent at = Ep^ a situation which confiicts with the standard WKB 
assumption that the potential can be approximated with a linear function /(A) = aX-\-b 
with a 7^ in the vicinity of the turning points. However, the quantum mechanical 
system is of course perfectly well defined independently of the validity of the WKB 
approximation. Let us now try to translate the instanton considerations of section 3.5.1 
to the quantum mechanical picture. These considerations of course refer to the situation 
where g > gc- By comparison of (3.172) and (3.195) we see that 

/I \ 2 f-U^+y) I 

AS(x = -^{x + y)^ = -j^ ^Vj,{X) - Ef dX, g>gc (3.196) 

and hence the matrix model probability that one eigenvalues tunnels through the bar- 
rier provided by the effective potential and escapes to infinity looks like the quantum 
mechanical probability that the fermion with energy Ep tunnels from the large po- 
tential well to the smaller one. Strictly speaking there is of course no probability for 
tunneling when the Fermi energy exactly coincides with the value of the potential at 
the local minimum. But this situation only holds when N is strictly infinite. For N 
not infinite we should replace the left hand side of (3.192) by 1 -|- It is easy to 
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show that this prescription causes the Fermi energy corresponding to the large well to 
be lifted by an amount /\Ef ~ O (y^^ which makes it agree exactly with the lowest 
level in the small potential well. Furthermore it appears that the level spacing is the 
same in the two wells. We thus have a situation which is remarkably similar to that 
of a symmetric double well. As for this case the degenerate energy levels will split 
into two and the lowest of these will have to be identified with the Fermi level. This 
corrected Fermi level, lies above the bottom of the small potential will and hence 

the fermion with energy E'p'^ can tunnel between the two wells. The probability for 
this process is e~'" where F is given to leading order in by (3.196). We see that in 
accordance with the stabilizing effect of the stochastic scheme there is no possibility 
for the fermion to escape to infinity. Furthermore we note that as was the case in the 
matrix model picture the tunneling is exponentially suppressed in the large- limit. 
Likewise the correction to the Fermi energy — Ep is exponentially small in N . 

So far we have restricted ourselves to the simple matrix model given by the poten- 
tial (3.168). As mentioned earlier for matrix models with potentials of degree larger 
than three the quantum mechanical system which arises by the stabilization procedure 
consists of A^ interacting fermions. However, in the limit N ^ oo these interacting 
systems are in fact equivalent to certain non interacting systems. A Hartree-Fock 
like approximation consistent with the large-A^ limit can be used to decouple the 
fermions [94]. With this Hartree-Fock like approximation the Fokker-Planck Hamil- 
tonian takes the same form as in (3.189) and (3.190). Of course the precise form of 
the Fokker-Planck potential Vfp depends on the original matrix model potential. For a 
matrix model potential of degree k the Fokker-Planck potential will have degree k — 2. 
It turns out that with the Hartree-Fock approximation the relations (3.193) and (3.195) 
still hold (in the region of the coupling constant space where the large- A^ limit of the 
original matrix model is well defined). Hence we get a quantum mechanical descrip- 
tion of all the matrix model m'th multi-critical points. We see that in the quantum 
mechanical picture the m'th multi-critical point is characterized by the (m — 1) first 
derivatives of the Fokker-Planck potential vanishing at (at least) one of the turning 
points corresponding to the Fermi energy. We furthermore get a quantum mechanical 
description of the criterion for stability of the original 1-matrix model. As mentioned in 
section 3.5.1 a m'th multi-critical model is unstable if we when we approach the critical 
point encounter a real instanton, i.e. a real zero for Wq*^"'(A) outside the support of the 
eigenvalue distribution. This situation in the quantum mechanical picture corresponds 
to the Fokker-Planck potential having a local minimum touching the Fermi level (for 
A^ = oo). 

3.5.4 The Double Scaling Limit 

Let us consider the situation where g > Qc (fig- 2a) and let us inspired by the calcula- 
tions on page 46-47 introduce a scaled variable z by 

/ _ 

A-x, = xJ^-^j z. (3.197) 
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Expressed in terms of z the Fokker-Planck Hamiltonian (3.190) reads 



^/P = "^/p(0; 5 - 5c) + x""}— — — h^^{z) 



3/2 

9c 



hfp{z) 



where Vjp(0; (/ — gc) is a second order polynomial in g — gc and 
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(3.198) 
(3.199) 

(3.200) 



(3.201; 



^9'J' 
NHg-g^y/^- 

We see that our problem of studying the Fokker-Planck Hamiltonian Hpp given by (3.189)- 
(3.191) has now been reduced to studying the scaled Hamiltonian hfp[z). For hfp[z) 
the double scaling parameter plays the role of i.e. the expansion parameter of the 
WKB approximation. Hence the double scaling expansion of our original matrix model 
is nothing but the WKB expansion of the quantum mechanical system given by hfp[z). 
We note that Vfp[z; \/ g — gc) contains terms which are subdominant in the double scal- 
ing limit. This is not unexpected since we know that a lot of terms which appear in 
the full expression for the free energy and the correlators of the matrix model vanish 
in the double scaling limit (cf. to section 3.1.5 and 3.2.3). Likewise it comes as no sur- 
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prise that we have in Hjp a trivial analytic term and that a factor of 
be extracted from the non analytic term in order to get a non trivial result in the 
double scaling limit. It follows from the conventional matrix model analysis bearing 
in mind that the Fermi energy of Hjp equals (^^ Tr ^) (cf. to page 51). The potential 
^fpi^i V9 ~ 9c) is of course of the same structure as Vfp. It has a global minimum as 
well as a local minimum. In the double scaling limit — )• (/c, — > 00, /i fixed the local 
minimum coincides with the local minimum of the potential 
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(3.202) 



The local minimum of Vds{z) is located at 2; = 1 and Vds{z = 1) = —2. In the WKB 
approximation the Fermi energy of the system given by hfp[z) is determined by 
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\I^S- Vfp{z)dz 



(3.203) 



where as usual zi[ef) and Zr{ef) are the left and right turning points for a classical 
particle moving in the potential Vfp[z) and where we have assumed that >> 1. The 
Fermi energy ej is related to Ep by 



Ef = Vfp{0;g - gc) + x'^ 



9 - 9c 



3/2 



(3.204) 
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Since the structure of the potential for the quantum mechanical system is unchanged 
by the scaling procedure it is obvious that ej coincides with the value of vjp at its local 
minimum. It can also be verified by direct computation. In particular one finds that 
in the limit N — )• oo, g ^ gc^fi fixed, ej = —2, i.e. the fermi energy equals the value 
of Vds at its local minimum. We note that Vds{z) also takes the value —2 for z = —2. 
We can carry over the symmetric double well analogy from the previous section to the 
scaled Hamiltonian as well. As before we should to be precise replace the N entering 
equation (3.203) by {N + |) and as before this causes the Fermi energy corresponding 
to the large well of vjp to be lifted by an amount which makes it agree exactly with the 
lowest level in the small well. Furthermore one finds that the level spacing is the same 
in the two wells. Again the degenerate energy levels will split into two and the actual 
Fermi energy 6^°'' will lie above the bottom of the small potential well. The probability 
that the fermion with energy tunnels between the two wells is now given by e~^"^'°-'- 
where 

2 relief) I 

^scai = T dzJvfp{z) - ef. (3.205) 

Here Zr{eY^) is the right turning point in the large well for a classical particle with en- 
ergy e^^ moving in the potential vjp and zi[eY^) the corresponding left turning point 
in the small well. Since the scaling prescription for A, (3.197), is nothing but the 
usual double scaling prescription (cf. to equation (3.178)) the Tscai appearing in equa- 
tion (3.205) is nothing but the double scaled version of the F appearing in section 3.5.3, 
i.e. the double scaled version of the instanton action, AS 

T,,ai = {AShs.i. (3.206) 

We hence have arrived at a quantum mechanical interpretation of the non perturbative 
ambiguity of the Painleve equation. In the matrix model picture we saw that the 
probability that the largest eigenvalue would tunnel through the barrier provided by 
the effective potential was suppressed in the large- limit but became of importance 
in the double scaling limit. The situation is very similar in the quantum mechanical 
picture. The probability that the most energetic fermion tunnels between the two 
potential wells of the Fokker-Planck potential is suppressed in the large- A^ limit but 
becomes of importance in the double scaling limit. Needles to say that the tunneling 
probability is a non perturbative quantity. 

Let us return for a moment to the expression (3.199) for the Hamiltonian of the 
quantum mechanical system. If we take the double scaling limit on operator level we 
are left with the following Hamiltonian 

hd, = -n^— + z^ -3z. (3.207) 

It is thus tempting to state that 2D quantum gravity in the double scaling limit is 
equivalent to a fermionic system with Hamiltonian hds- However, some clarification is 
needed at this point. The operator h^s is not automatically self-adjoint. To ensure self- 
adjointness its domain must be carefully specified. As explained in references [97, 98] 
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one can easily by means of the WKB approximation convince oneself that there exists 
a one parameter family of self-adjoint versions of h^s all of which have discrete spectra. 
Actually it can be proven (also without referring to the WKB approximation) that 
this is the situation for a large class of Hamiltonians with potentials unbounded from 
below, namely Hamiltonians, /i, of the type 

h = -—+v{y) (3.208) 

where v[y) — )• — oo for y — )• — oo so fast that 

dy 



< oo (3.209) 

^{y)\ 

and where v[y) behaves in the same way when j/ — )• oo or v[y) — )• oo for y — )• oo. We 
refer to the mathematical literature, for instance [99] for a rigorous discussion of these 
points. Let us just note that the condition (3.209) implies that a classical particle, if it 
is not trapped in a local minimum of v[y) will move to — oo or possibly +00 in a finite 
time. Furthermore let us note that the condition (3.209) hints to us the quantization 
of energies since in regions where the WKB approximation is valid we have for the 
density of states 

p{e)= I ^^^^e{e-v{y))dy. (3.210) 

■' 2TrJe-v{y) 



As shown in reference [98] in the energy region below but close to the Fermi level of 
the full Vfp the level spacing is the same for all the different self-adjoint extensions of 
hds- However, any two extensions will in this region differ by an overall displacement of 
the energy levels. This allows us to advocate exactly one of the self-adjoint extensions 
of hds as the one relevant for 2D quantum gravity, namely the one having an energy 
level coinciding with the Fermi energy of the full Fokker-Planck potential, vjp. Only 
with this Hamiltonian we will reproduce correctly the non perturbative effects of our 
original model. 

Fven though the quantum mechanical system given by (3.207) is hence perfectly well 
defined the calculations we were originally aiming at being able to carry out, namely 
calculations of expectation values of operators like Tr can not be performed 
without referring to the full Fokker-Planck potential. This is of course closely related 
to the fact that for these expectation values a multiplicative renormalization by some 
power of {g — gc)~^^^ was needed before a double scaling expansion could be extracted 
(cf. to section 3.2.2). With the Hamiltonian (3.207) we can only access quantities which 
can be expressed entirely in terms of the double scaling parameter. These aspects are 
discussed in more detail in reference [98] where it is also shown how one can easily 
extract the correct scaling behaviour of observables from h^s just by introducing an 
appropriate cutoff in the negative z. 
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3.5.5 A Strong Coupling Expansion ? 

The considerations outlined in the previous section seem to hint to us the possibihty 
of a convergent strong couphng expansion of 2D quantum gravity. To explain why let 
us recall a well understood problem from quantum mechanics, namely the anharmonic 
oscillator 

H = -—- + x^+gx\ g>0. (3.211) 

For the operator H the perturbative expansion in powers of g is not a convergent 
power series expansion but only an asymptotic expansion. The term gx'^ is not a small 
perturbation for any g > and clearly g — )• —g changes drastically the nature of H. 
However, the anharmonic oscillator has a convergent strong coupling expansion. This 
can be seen by a simple scaling argument. Setting 

y = 9'^'x (3.212) 

we find 

H{x) = \-^'^H{y) (3.213) 

where 

H{y) = -'^,^y'^Xy\ A = ^. (3.214) 

There exists a theorem stating that a term v{y) is an analytic perturbation of an 
operator Ho{y) = + Vo{y) if 

1. D{v) D D{Ho) 

2- I bV'l I2 ^ «| l-f^oV'l I2 + ^1 IV'I I2 for some a^b ^ R and all t/j G D[Ho). 
The potential v[y) = y^ clearly satisfies the criteria above for being an analytic per- 



turbation of Ho{y) = —j^ + y*- For an eigenvalue En of the original operator H we 



can therefore write 

00 

EMX)) = X-'''J2^nkX'' (3.215) 

k=0 

where the power series has a finite radius of convergence — and similarly for the 
eigenfunctions. The situation is very similar in the case of the Hamiltonian hfp[z). 
The expansion in powers of is only an asymptotic expansion. However, we can by 
means of a simple scaling prescription transform the Hamiltonian hfp[z) into one for 
which the natural expansion parameter involves negative powers of h. If we introduce 
scaled variables by 

z = n^/^y^ X = 4j, = — (^Y' (3-216) 

we can write 

hfpiz) = X-^l^hfp{y), (3.217) 
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hfpiy) 



J2 



(3.218) 



where 



Vfp{y,X,N) = y^-3Xy + 



2 

2x2/5 

iV 



3A 
3A 



?/ + -y 

2 4 

— V H — V 

2 4 



(3.219) 



and the double scaling limit is obtained as before by letting — )• oo, — )• (/c, but A 
fixed. The non-trivial information is contained in the Hamilton function hfp{y) which 
in the double scaling limit goes to 



hds{y) 



dy"^ 



+ Vds 



dy"^ 



+ y^- 3Aj/ 



(3.220) 



and we have formally the same situation as for the anharmonic oscillator: — 3Aj/ looks 
like a small perturbation with respect to y^. It is thus natural to expect a strong 
coupling expansion of the energy eigenvalues of hds{z) of the form 



k=0 



k=0 



(3.221) 



where the series is convergent — and similarly for the eigenfunctions. The domain of 
the perturbation v[y) = — 3Aj/ does not include the domain of ho[y) = —d^/dy^ + y^ 
so we have no rigorous proof of this conjecture, but one should keep in mind that the 
requirements (1) and (2) mentioned above are only sufficient conditions, not necessary 
conditions, for analyticity. 
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4 D>2 



4.1 Analytical Approaches 

What caused the progress in our understanding of 2-dimensional quantum gravity was 
the discovery that the matrix model (3.1) exactly generates 2D orientable simplicial 
manifolds weighted with the Einstein Hilbert action. It is natural to look for a simi- 
lar mapping of the partition function of 3D and 4D simplicial gravity given by (2.11) 
and (2.13) onto an exactly solvable model. In two dimensions one automatically gets 
a simplicial manifold if one glues together J-simplexes along their [d — l)-dimensional 
sub-simplexes until one has a closed complex. This is not the case in 3 and 4 dimen- 
sions. Constructing a 3- or 4-dimensional simplicial complex, T, following the recipe 
just given one only obtains a so-called pseudo-manifold. A pseudo-manifold, T, is a 
simplicial manifold if and only if x{T) = 0. In general x{T) > 0. This complica- 
tion implies that one needs more parameters than the two required by the Einstein 
Hilbert action in order to construct an acceptable model of simplicial gravity in 3 and 
4 dimensions. One way to search for viable theories of 3D and 4D quantum gravity 
is by trying to generalize the construction described on page 11 [45]. The most naive 
generalization of the construction does not work since only two parameters appear. 
There exist other generalizations involving more parameters [100, 101, 102]. However, 
none of these models are capable of at the same time eliminating complexes that are 
not simplicial manifolds and restricting the topology, let alone providing us with a 
topological expansion. Lacking a topological classification of 3D and 4D manifolds 
we do of course not really know what kind of expansion possibilities we can expect 
to encounter, if any. It is not an unreasonable line of action to try to generalize the 
two-dimensional models. However, since 3- and 4-dimensional geometry is much more 
complicated than two-dimensional geometry it might be necessary to consider models 
involving more complex structures. 

A model of three-dimensional simplicial gravity involving the structure of quantum 
groups has been suggested by D. Boulatov [103, 104]. To define the model one uses 
as the starting point a real function, </, of three variables defined on the quantum 
group SUq{2) where q = _ The dynamical variables of the model are the 

fourier coefficients of the function g. The free energy of the model can be written as 
a sum over all connected orientable 3-dimensional pseudo-manifolds. A given pseudo- 
manifold, r, appears with a weight which in addition to the Einstein Hilbert weight 
factor involves the factor /g(r)A^(-^) where 



and where Iq{T) is a topological invariant of T, the so-called Turaev Viro invariant 
associated with the quantum group SUq{2) [105]. If T is a triangulation of an orientable 
3-dimensional manifold, Af, it holds that [106] 






I,{T) < AJW 
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where h[M) is the Heegaard genus of the manifold M. The Heegaard genus, h{M), 
equals zero if and only if M is homeomorphic to S"^ . Otherwise h[M) > 1. For the 
three-sphere one has Iq{S^) = 1. Hence if one could send Ag to zero keeping valid 
the inequality (4.223) one would be able to obtain a situation where only simplicial 
manifolds homeomorphic to would contribute to the free energy of the model. How- 
ever, such a limiting procedure does not immediately make sense, since Ag — )• implies 
q ^ and q is supposed to be a root of unity. 

Ooguri has suggested another line of action for how to obtain from Boulatovs model 
the free energy of 3-dimensional simplicial Einstein Hilbert gravity where only simplicial 
manifolds homeomorphic to the sphere are taken into account [107]. This line of action 
is very similar to the line of action followed in section 3.1.5 where the free energy of 
the hermitian 1-matrix model was calculated. There we first calculated the 1-loop 
correlator using the Dyson Schwinger equations and afterwards integrated the result 
to obtain the free energy. Similarly Ooguri suggests that one should go via the Dyson 
Schwinger equations to obtain the free energy of 3-dimensional simplicial quantum 
gravity. Ooguri defines a 3-dimensional analogue of the matrix model expectation 
value ( Tr (f>^) and derives a Dyson Schwinger equation analogous to equation (3.14) 
for this quantity. Just as the expectation value ( Tr (f>^) in the matrix model case can be 
written as a sum over 2-dimensional simplicial manifolds with a loop of length n inserted 
Ooguris expectation value can be written as a sum over 3-dimensional pseudo-manifolds 
with a 2-dimensional surface of a certain area and a certain number of handles inserted. 
And as in the matrix model case one can derive the factorization condition (3.16) which 
makes it possible to arrange that only 2-dimensional manifolds of genus zero enter the 
Dyson Schwinger equation Ooguri can write down a set of factorization rules by means 
of which one can arrange that the Dyson Schwinger equations of the 3-dimensional 
model involve only 3-dimensional manifolds homeomorphic to S^. For this version of 
the Dyson Schwinger equations the solution related to the insertion of the 2-dimensional 
surface constituted by the boundary of a 3-simplex is a simple derivative of the free 
energy of 3-dimensional Finstein Hilbert gravity for simplicial manifolds homeomorphic 
to S^. In principle it should be possible to calculate explicitly the observables of this 
model by means of the Dyson Schwinger equations and the factorization rules. However, 
such calculations have not yet been carried out. In the next section we will describe 
an approach to 3-dimensional simplicial quantum gravity where observables have been 
calculated, namely numerical simulations. Let us close this section by mentioning that 
Boulatovs model has also been generalized to 4 dimensions [108]. 

4.2 Numerical Simulations based on DT 
4.2.1 The Numerical Method 

As appears from the previous section so far no powerful analytical techniques for study- 
ing the partition function (2.18) of dynamically triangulated gravity for d > 2 have 
been developed. However, the 3- and 4-dimensional models can be studied by numer- 
ical means. The strategy when studying these models numerically consists in using 
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Monte Carlo techniques to generate J- dimensional simplicial manifolds homeomorphic 
to S'^ with weights given by the Einstein Hilbert action. For the generated simplicial 
manifolds one measures certain quantities and averaging over a large number of mea- 
surements one obtains the quantum averages of the quantities considered. To generate 
the simplicial manifolds one starts out with a given triangulation of S'^ and moves 
around in the space of combinatorially equivalent triangulations with a set of so-called 
moves, i.e. local changes of the triangulation^. These moves must fulfill the requirement 
of ergodicity, i.e. they must allow us to reach any simplicial manifold combinatorially 
equivalent to the initial one. Furthermore they should of course not change the topol- 
ogy of the manifold. In the grand canonical ensemble, i.e. the ensemble of manifolds 
where the volume is allowed to vary an ergodic set of moves has been known since 
1930. These are the so-called Alexander moves [109]. The Alexander moves are, how- 
ever, not convenient from a numerical point of view. More convenient are the so-called 
(p^q) moves. In d dimensions there are [d -\- 1) such moves. They can be described in 
the following way. Let us consider a [d -\- 1) simplex, n^+i. Its boundary, drid+i con- 
sists of {d-\- 2) J-simplexes. The recipe for performing a (p, q) move in a J-dimensional 
simplicial manifold reads: Pick out a collection of p J-simplexes which can be mapped 
onto a subset S of drid+i and replace it by the complement S in drid+i consisting of 
[d-\-2 — q) J-simplexes. The (p, q) moves were shown to be equivalent to the Alexander 
moves for d = 3 in reference [17] and for J = 4 in reference [110]. As opposed to what is 
the case in two dimensions in 3 and 4 dimensions no set of moves which are ergodic in 
the canonical ensemble, i.e. the ensemble of manifolds with a fixed volume, is known. 
One is hence forced to perform simulations in the grand canonical ensemble. However, 
there exists a method by means of which one can imitate a fixed volume simulation 
while still ensuring the ergodicity of ones algorithm [HI]. As explained in section 2.1 
a necessary condition for the partition function given by (2.18) to make sense is that 
the number of triangulations with the topology of S'^ and a given volume, N^, must be 
exponentially bounded by the volume. Let us for a moment assume that such a bound 
exists, i.e that the partition function behaves as in (2.25). Normalizing the partition 
function such that the entropy factor for manifolds consisting of J-simplexes equals 
one we can write 

Z{k,, k,_2) ~ E e-('^<^-'^5('^<^-)^'')(^<^-^S) (4.224) 

where 

>^U>^d-2r'' = «^(«.-2) + (4.225) 

We see that fixing at the value K'^^{Kd-2Y^ ^ corresponding to some not too small 
volume and adding a term /\nd\Nd — A^°|, < Ak^ << 1 to our original action 
results in the volume of the manifolds generated being peaked at and distributed 
around this value according to 

P{Nd) ~ exp(-A/c,|iV, - iV°|). (4.226) 



®Two triangulations are combinatorially equivalent it they can be subdivided into the same finer 
triangulation. 
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This observation is very useful. First it makes it possible to check whether the expo- 
nential bound referred to above exists or not and if so to investigate further the nature 
of its subleading correction. Let us briefly describe how this works in practise. To 
check whether an exponential bound exists or not we start by choosing a given volume, 
A^°, a value of K,d-2 and a value of Ak^- Then we search for a value of k.^ for which 
the actual volume of the generated manifolds fluctuates in an interval around N^. If 
one succeeds in finding a value of k.^ for which such a situation is encountered one has 
a first estimate of ^"^^{^01-2)^^ ^ ■ By examining the distribution of Nd around one 
can improve the estimate for K'^^{Kd-2Y^ ^ ■ (If i^d does not equal K'^^{Kd-2Y^ ^ the slope 
of In [P[Nd)) will be different above and below and this difference in slope tells us 
by which amount we should correct our initial guess.) The condition that the entropy 
of manifolds is exponentially bounded reads 



By repeating the fine-tuning process described above for different values of one can 
test if this condition is fulfilled. In particular one can investigate the nature of the 
(hopefully) subleading correction to the exponential bound. We note that if for d > 2 
we have a ^/str like in 2 dimensions i.e. f{Nd,Hd-2) = N] ~^ the equation (4.225) 
reads 



and a measurement of the distribution of the volume of our manifolds would hence 
allow us to determine 7(/€d-2)- 

Fixing to its critical value ^"^^{^01-2)^^^ and adding to the action the term /\Kd\Nd— 
N^\ is also what allows us to imitate a fixed volume simulation without violating the 
ergodicity of our algorithm. By adjusting Ak^ we can obtain an appropriate width 
of the distribution of volumes and by making measurements only on manifolds with 
volume we effectively study the fixed volume partition function. The value of Ak^ 
is typically only a few percent of K'^^{Kd-2Y^ ^ ■ The strategy of fixing Hd at its critical 
value is fully compatible with our desire to search for a critical point where a contin- 
uum theory might exist since as explained in section 2.1 a necessary ingredient in our 
recipe for defining a continuum limit is to let Hd approach from above. Let us point 
out again that the aim of the numerical simulations is to search the restricted coupling 
constant space Hd = '«d('«(i-2, {gi]) for a critical point where an appropriately defined 
susceptibility for quantum gravity diverges. Here the coupling constants {gi] might 
be coupling constants for matter fields or for higher derivative terms. We note that in 
order to pursue this aim one has to perform a fine tuning of Hd for each choice of the 
remaining coupling constants. In addition we have chosen to perform for each choice 
of coupling constants yet another fine tuning, namely we introduce a probability for 
performing each of the [d moves which we adjust so that in a given simulation all 
the [d -\- 1) different moves are performed approximately the same number of times. 
The purpose of this arrangement is to avoid bias towards special types of manifolds 




00. 



(4.227) 





62 



caused by a difference in acceptance rates for tlie different types of moves. Let us aiso 
mention tliat we restrict tlie o(ni)'s in tlie foiiowing way 

o{n,)>d-i + l, i<d- 2. (4.229) 

Of course we liave 0(^^-1) =2. To do tlie updating of tlie geometry we use a standard 
Metropolis algorithm. Possible matter fields living on the geometries are updated 
using well established techniques from fixed lattice simulations. However, special care 
must of course be taken for the updating of matter fields related to changes in the 
geometry. We will not enter into a discussion of these points here. A description of 
the line of action taken in the simulations that we will report on below can be found 
in references [112] and [113]. 



4.2.2 Observables 

In quantum gravity we do not have the same type of observables as in conventional 
quantum field theory. Expectation values of operators which refer to specific points 
in the space time manifold have no meaning. Observables which do make sense are 
averages of integrated local operators. One such observable which we have made use 
of in our simulations is the average integrated scalar curvature per volume 



(R) = 7 .,rr- ^ ^ - Kd+1,2 (4.230) 



d 



where the notation is as in section 2.1 (cf. to equation (2.7) and (2.10)). Another type 
of observable which is available in quantum gravity is the average correlation between 
local operators 0{^i) and 0(^2) defined as a function of geodesic distance, for example 
as 

{G{r)) = {[ [ ^/'6A/^KeO^'6AA(6)0(6)0(6)<^(<6,6)-r)) (4.231) 



M JM 



where (^(1^1,(^2) is the geodesic distance between i^i and ^2 for a given Riemannian 
manifold. A particularly simple example of such an observable is the average volume 
{V{L)) of the ball with radius L defined as the volume of the set of points {^i} at 
geodesic distance d[^i^^o) < L from some reference point, i^o- By means of {V{L)) we 
can define the Hausdorff dimension, J//, describing the fractal structure of space time 
by 

{V{L)) L^i^. (4.232) 

In numerical simulations of simplicial gravity we obviously have the possibility of study- 
ing in detail the fractal structure of our universes. From a numerical point of view it 
is convenient in stead of working with the strict definition of geodesic distance on a 
simplicial manifold to work with the so-called n^-distance between J-simplexes. We 
define the n^-distance, D, between two J-simplexes in a J-dimensional simplicial mani- 
fold as the length of the shortest path between the two J-simplexes obtained by moving 
between centers of neighbouring J-simplexes. By (D) we mean the average value of 



63 



D for a given manifold averaged over our ensemble of manifolds. Hence (D) can be 
thought of as a typical radius of our universes. We have measured the average value 
of D as well as its (average) distribution, P{D). From these measurements we can 
in principle extract du by fitting P{D) to the functional behaviour P{D) ~ D'^^~^. 
Likewise we can in principle determine the so-called cosmological Hausdorff dimension 
dch defined by [45] 



{D) ~ N,_ 



(4.233) 



In practise it is hard to extract in a reliable way any of the Hausdorff dimensions. (We 
note that the two definitions are not necessarily equivalent.) However, (D) as well as 
(R) have proven very convenient when scanning the coupling constant space of the 
models of dynamically triangulated gravity. 

To investigate further the nature of a possible phase transition we will need a 
susceptibility like observable. One candidate for such an observable is the curvature- 
curvature correlation function. In continuum notation it reads 

x{K,_,) = {j j d'i,d'i,^Jg{^)R{i,)^Jg{^)R{i,))-{j d'i^)R{()f (4.234) 

and it is easily seen that it is proportional to the second derivative of — log Z with 
respect to the inverse gravitational constant (cf. to equation (2.4)). In the context of 
dynamically triangulated gravity with fixed x{>^d-2) turns into 



xiKd-2,m) = {Ni,)N,-{m-2) 



(4.235) 



What we hope to encounter is a point where x{^^d-2^ Nd) I Nd diverges when Nd — ^ oo. 
We have in addition to x{^d-2^ Nd) also studied another somewhat related quantity, 
namely (i?^) — {RY where 



(7^ E <^d- 



Cd - o{nd-2) 
o{nd-2) 



Kd+1,2). (4.236) 



This expression for the average squared curvature of a J-dimensional simplicial manifold 
comes about when one assigns to each of its ( J — 2)-simplexes a volume density Ki^_2 = 



nd-2 



{cd - o{nd-2))IVr. 



nd-2 ■ 



We will di 



iscuss 



o{nd-2) I Kd+1,2 and a curvature density R 
this construction in more detail in section 4.4.3. The correlator (i?^) — {RY has turned 
out to be a useful indicator for changes in the geometrical system. 

In the case where matter is coupled to gravity one of course has additional observ- 
ables associated with the matter fields. As for the geometry possible observables for the 
matter systems are restricted to global quantities and correlation functions defined in 
terms of geodesic distances. In the following we will only discuss aspects related to the 
behaviour of global observables of the matter systems. For a discussion of the problem 
of defining correlation functions as well as the problem of extracting from these critical 
indices we refer to [114]. 
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4.3 Numerical Results for D=3 



4.3.1 Pure Gravity 

The justification for simuiating 3-dimensionai dynamicaiiy trianguiated gravity is mainiy 
tliat it is a naturai intermediate step between tlie exactiy soivabie 2-dimensionai case 
and tlie far more complicated 4-dimensional case. Classical 3-dimensional general rel- 
ativity it trivial. However, the Einstein Hilbert action is in three as in four dimensions 
unbounded from below and non-renormalizable. Hence it might very well be that nu- 
merical simulations will reveal to us a theory with a non trivial structure. In particular 
coupling to matter might lead to complex interactions. 

The first step when aiming at simulating three dimensional simplicial quantum 
gravity is of course to check that the partition function of the model makes sense. In 
reference [16] it was verified by numerical means that the number of triangulations of 
the three-sphere with a given volume is exponentially bounded by the volume. Having 
settled this question one knows that the line K3 = ^^(ki) described in section 2.1 
and 4.2.1 exists and one can start looking for critical behaviour along this line. Several 
groups have pursued this idea [17, 115, 116, 117] and their results agree. 

One finds that when moving along the line K3 = ^^(ki) one encounters at a certain 
value of Ki (between 3.5 and 4.0) a seemingly sudden change between two types of 
geometry. The two regions in the coupling constant space corresponding to the two 
types of geometry are denoted as the hot and the cold phase of 3-dimensional simplicial 
quantum gravity. For small values of Hi the geometrical system is in the hot phase. The 
hot phase is characterized by the space time manifolds being very crumpled. Typically 
there will be vertices of very high order and the number of vertices increases only very 
slowly with the volume. The average integrated curvature per volume (R) will be small 
and eventually for small enough values of Hi become negative. The average radius of 
the universes (D) will also be small and almost volume independent. In addition the 
distribution of D, P{D) will be strongly peaked. The Hausdorff dimension defined 
in either way will hence be very large. Possibly it is even infinite. Presumably there 
is no 7str in this phase [16, 17]. The cold phase comes into existence when Hi gets 
large. In this phase the space time manifolds are very extended. Their average radius, 
(D), is large and grows almost linearly with the volume. Furthermore there are large 
fiuctuations in (D) and P{D) is a very broad distribution with an almost plateau like 
maximum. The Hausdorff dimension of the space time manifolds in the cold phase is 
probably close to one. The average curvature (i?) is large and positive. It seems that 
in this phase a ^str does exist [17]. The cold phase has been interpreted as signaling the 
dominance of the conformal mode. Being regulated the Fuclidean path integral can of 
course not diverge. However, the 1-dimensional space time manifolds with a rapidly 
changing connectivity are probably the lattice analogues of the conformal instability. 
For small values of Hi the conformal instability is suppressed by the large entropy of 
very compact universes. In connection with the transition between the hot and the 
cold phase pronounced hysteresis is observed. This is illustrated in figure 3 where we 
have shown the behaviour of (i?) as our geometrical system goes through a thermal 
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cycle. 




Figure 3. Hysteresis curves for 3- dimensional dynamically triangulated 
gravity. Triangles: pure gravity. Circles: gravity plus one Ising spin 
system when the coupling between matter and geometry is maximal, i.e. 
/3 = Ns = 10000. From reference [119]. 

The fact that hysteresis is observed shows that the transition between the two phases 
is of first order. Hence we have no hope of defining along the lines of conventional 
statistical mechanics a continuum theory of gravity with a propagating massless par- 
ticle. It seems as if the classical result that no gravitational waves can propagate in 3 
dimensions survives quantization. The situation might be different if matter is coupled 
to gravity. We will consider this situation in the next section. Let us mention here 
that in the static triangulation approach to simplicial quantum gravity one finds for 
d = 3 a second order phase transition both in the case of pure gravity as well as in 
the case where a term is added to the Finstein Hilbert action [118]. So far no 
clear understanding of this rather fundamental difference between the two models of 
quantum gravity has been achieved. 

4.3.2 Gravity plus Matter 

Studying the interaction between space time and matter is of course interesting in its 
own right. For d = 3 we furthermore have the possibility that the presence of matter 
fields might cause the first order phase transition observed for pure gravity to change 
into a (from continuum physics point of view) more interesting second order phase 
transition. 
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So far only the case of 3-dimensional simplicial gravity coupled to an Ising spin 
system has been investigated. As mentioned in section (2.2), when coupling matter 
to simplicial gravity one can choose to place the matter field variables either on the 
vertices of the J-dimensional simplicial manifold or in the center of its J-simplexes. We 
have chosen the latter possibility [119]. In this way we are sure that the number of 
matter field variables grows with the volume of the universe. This is not necessarily 
true if the matter fields variables are placed on the vertices of the simplicial manifold 
(cf. to section 4.3.1). Our partition function reads 

Z{Ks,K,,(3) = j:e-^^^ J: 1 e-^^^e'^^(..)(^--^) (4.237) 

Ns T~53(Af3) ^ 'v-' ) {a} 

where (Ti G { — 1,+1} and where J2(i,j) usual denotes the sum over neighbouring 
pairs of spins. We note that when the spins are placed in the center of the 3-simplexes 
each spin has 4 neighbours. By J2{a} we mean the sum over all spin configurations for 
a given manifold. As usual we are only interested in exploring the restricted coupling 
constant space K3 = k,'^[k,i^ jS). Hence we have effectively only two coupling constants 
Ki and /3. What we aim at is mapping out the phase structure of the model in the 
(/€i,/3) plane. To prepare for the presentation of the result it is convenient first to 
recall the characteristics of the Ising model on a regular lattice. A one dimensional 
Ising chain never becomes critical. The spins are always disordered. For d = 2 the 
Ising model has a second order phase transition from a disordered to an ordered phase 
at some value /3c of /3 and mean field theory indicates that a similar behaviour is seen 
for d > 2. 

The phase diagram of the model (4.237) is shown in figure 4. The dotted line 
is the line Hi = Ki(/3), i.e. it indicates the transition between the two phases of the 
geometrical system. (Due to hysteresis this line is not unambiguously defined. We 
shall explain shortly what we mean by /€j(/3).) The almost horizontal line indicates 
a transition between a non magnetized and a magnetized phase for the Ising model. 
The location of the line /3 = /3c('«i) is determined with the aid of the magnetization 
curves for the spin system (for N3 = 10000). The transition between the two phases 
becomes sharper as the volume of the universe increases and we presumably face here 
a second order phase transition. We note that the dependence of /3c on Hi is rather 
weak. In the elongated phase of the geometrical system we see no sign of a phase 
transition for the spin system. Presumably in the infinite volume limit the spin system 
has zero magnetization for all values of /3. (We refer to reference [119] for details.) The 
absence of a phase transition for the spin system in the cold phase of gravity seems 
to show that the Hausdorff dimension of our space time manifolds is the dimension 
relevant for coupling to matter. Accordingly it is tempting to conjecture that for the 
phase transition of the spin system encountered in the hot phase of gravity, where the 
Hausdorff dimension is very large, the critical exponents take their mean field values. It 
would be interesting to test this conjecture but a reliable determination of the critical 
exponents of the matter system requires more extensive numerical studies. Figure 4 
shows that there is a clear back reaction from the spins on the geometry. The spins 
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affect the geometry by shifting k.^ towards smaiier vaiues. The effect is iargest when the 
spin system is itseif criticai and talces its pure gravity vaiue for /3 — )• and /3 — )• oo. 
This is what one wouid intuitiveiy expect. When /3 — )• or /3 — )• oo the Ising spins are 
independent and individuaiiy they can not inffuence the geometry. Oniy when there 
are correiations between the spin variabies their presence might be of importance. It is 
interesting to note that the spin system seems to be counteracting its own criticality. It 
effectively pushes the geometry towards a lower Hausdorff dimension. A similar effect 
is seen in 2 dimensions when multiple spins are coupled to gravity [37]. 
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3.9 4.2 

Figure 4: The phase diagram for 3-dimensional dynamically triangulated 
gravity coupled to one Ising spin system. Dotted line: the transition be- 
tween the hot and the cold phase of the geometrical system. Full drawn 
line: transition between a disordered and an ordered phase of the spin sys- 
tem. Filled circles are results for determined for N3 = 10000. From 
reference [119]. 

The back reaction from the spins does not change the order of the phase transition for 
the geometrical system, however. For all values of /3 one observes hysteresis in connec- 
tion with the transition between the hot and cold phase of the geometry. Although the 
hysteresis gets less pronounced as /3 — )• /3c it never completely disappears. In figure 3 
we have shown the behaviour of (R) as the geometrical system for /3 = /3c('«i) goes 
through a thermal cycle. When an Ising spin system was coupled to two-dimensional 
dynamically triangulated gravity one found that at the critical point of the spin system 
the exponent 70 describing the fractal structure of space time was changed from — | 
to — |. In the 3-dimensional case the spin system can only become critical when the 
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geometry IS in its hot phase and as nientioned earher m this phase presumably no ^str 
exists. Along the line Ki = kI{J3) we might have a ^fstr <iiid it would be interesting to 
explore what happens to this index as well as the critical indices of the spin system 
exactly at Hi = Needless to say that this is a project which calls for a consider- 

able amount of computer time. Let us finally comment on our definition of k.^ . If one 
starts out with values of Hi and /3 corresponding to the hot phase of the geometrical 
system and changes Hi in small steps keeping /3 fixed the value of (R) will trace out 
the lower branch of a hysteresis curve as in figure 3. For values of Hi well inside the 
hot phase the trajectory is unique and we define the critical value of Hi as the value 
of Ki for which this trajectory when extrapolated intersects the parts of the hysteresis 
curves corresponding to the cold phase of the geometrical system. The filled circles in 
figure 4 indicate values of k.^ obtained according to this philosophy using universes of 
size Ns = 10000. 

Since hysteresis occurs even at /3 = /3c coupling an Ising spin system to 3-dimensional 
gravity does not improve our possibilities of obtaining from the lattice approach a con- 
tinuum theory of quantum gravity in three dimensions. One could of course try to 
couple more spin fields or other types of matter fields to the geometrical system but 
this line of action has so far not been pursued. Let us mention here that 3-dimensional 
dynamically triangulated gravity coupled to an Ising spin system has also been in- 
vestigated in references [120] and [121]. In reference [120] the spins were placed on 
the vertices of the simplicial manifolds while in reference [121] both the case of spins 
living on the vertices and the case of spins living in the center of the J-simplexes were 
explored, but only for values of Hi corresponding to the hot phase of the geometrical 
system. In neither reference the possible infiuence of the spin system on the phase 
transition of the geometrical system was investigated. 

4.4 Numerical Results for D=4 
4.4.1 Pure Gravity 

Numerical investigations have shown that the number of triangulations of S'^ with a 
given volume is exponentially bounded by the volume. Hence the partition function of 
4-dimensional dynamically triangulated gravity is well defined in a certain region of the 
coupling constant space, K4 > k,1[k,2). As already explained, when aiming at defining a 
continuum version of the model one must search the line K4 = k,1[k,2) for critical points. 
Such investigations have been carried out by several groups [18, 19, 20, 21] and their 
results are in agreement. 

When one explores the restricted coupling constant space K4 = k,1[k,2) one finds that 
the geometrical system can exist in two phases. For small values of k,2 the geometrical 
system is in a hot phase and for large values of k,2 the system is in a cold phase. The hot 
and the cold phase of 4-dimensional simplicial gravity have the same characteristica 
as the hot and the cold phase of 3-dimensional simplicial gravity (cf. to section 4.3.1). 
However, the nature of the transition between the two phases is different in 3 and 4 
dimensions. In 4 dimensions no hysteresis is observed. This opens the possibility that 



69 



the phase transition might be of second order and hence that it might be possible 
to define in the vicinity of the critical point a continuum theory with macroscopic 
correlations. In figure 5 we have shown the behaviour of the susceptibility x('^2)/^4 as 
the system passes from one phase to the other. The susceptibility clearly has a peak 
and the position of this peak coincides with the value of k,2 for which the change in 
geometry is seen to take place. The height of the peak does increase with the volume 
of the universe. However, our data does not allow us to determine whether we face a 
divergence or not. In connection with the the transition between the hot and the cold 
phase we also observe a peak in the correlator (i?^) — (i?)^. Neither for this correlator 
have we been able to analyse in detail the volume dependence. It seems unlikely that 
the transition should be of first order but it is not excluded that it could be of order 
higher than two. 




Figure 5. Four-dimensional dynamically triangulated gravity. The be- 
haviour of the susceptibility x('^2)/^4 as the geometrical system undergoes 
the transition from the hot to the cold phase. N4 = 4000(«), N4 = 9000(n) 
and N4 = 16000(x). From reference [123]. 

Fven if we succeed in arriving at a convincing argument that the transition is of second 
order there remains a number of unclarified points. A natural part of a prescription 
for defining a continuum limit is to introduce a scaling parameter, a, to be considered 
as the lattice spacing, and to send a to zero and N4 to infinity while keeping fixed the 
physical volume of space time a'^N4. This implies the following scaling for (R) 

{Rcontinuum) {Rlattice) 7^ • (^4.238^ 



70 



Hence it seems that we are led to the conclusion that the curvature of our continuum 
space time manifolds can only be finite if the lattice curvature scales to zero at the 
critical point. However, one finds that the value of (Riattice) hes relatively far above 
zero at the critical point. (The same is true in 3 dimensions (cf. to figure 3).) It is 
still not completely excluded that (Riattice) 7^ is a finite volume effect but we 
consider it unlikely. We are convinced that (Riattice) 7^ is not due to the missing 
tessellation of fiat space for d > 2 since already when averaged over only the 
[d — 2)-simplexes belonging to a single J-simplex the curvature can easily take the 
value zero. One point of view to take is that the pure gravity theory is too naive and 
that (Riattice) ^iU actually be zero for a more realistic theory of gravity including 
for instance matter fields of higher derivative terms. Another point of view would be 
to give up the interpretation of (Rcontinuum) as the actual physical curvature of our 
continuum space time manifold. Lacking an alternative interpretation we have chosen 
to investigate a class of more elaborate models of dynamically triangulated gravity in 
4 dimensions. These investigations are reviewed in the following two sections. There 
is another observable whose value at the critical point has important consequences for 
a possible continuum limit, namely the Hausdorff dimension. As mentioned earlier 
in the hot phase of dynamically triangulated gravity the Hausdorff dimension is very 
large, possibly infinite, while in the cold phase it is close to one. Since the transition 
from large to small values of the Hausdorff dimension is continuous one might find a 
Hausdorff dimension close to four at the critical point. However, we will refrain from 
making any claims abut the value of dn or dch at the phase transition point since we 
do not believe that the data allows for a reliable estimate. 

Four-dimensional simplicial quantum gravity has been studied also in the static 
triangulation approach [122]. In this approach one finds a first order phase transition 
for the model describing pure gravity. However, the nature of the transition changes 
when a higher derivative term is added to the action. We will return to this point in 
section 4.4.3. 

Let us finally mention that in reference [13] the effect of choosing a measure different 
from the uniform one in dynamically triangulated 4-dimensional gravity was investi- 
gated. It was found that a change of the measure changes the location of the critical 
point. The abruptness of the transition between the hot and the cold phase was in gen- 
eral weakened when the measure was chosen different from uniform. The value of (R) 
at the critical point did not change significantly. A future aim of these investigations 
would be a detailed comparison of appropriately defined critical exponents. 

4.4.2 Gravity plus Matter 

Apart from the study of the interaction between matter and space time being interesting 
in its own right we also have the possibility that the introduction of matter fields in 4- 
dimensional dynamically triangulated gravity will improve the scaling of the curvature 
at the critical point (cf. to section 4.4.1). We have studied 4-dimensional simplicial 
quantum gravity coupled to an Ising spin system as well as 4-dimensional simplicial 
quantum gravity coupled to various numbers of Gaussian scalar fields. In all cases the 
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matter field variables were placed in the center of our J-simplexes. Here we will only 
briefly review the results. For details we refer to [ffS]. 

Our partition function for the system consisting of one Ising spin system coupled 
to 4-dimensional dynamically triangulated gravity is completely analogous to (4.237) 
and the phase diagram of the model looks precisely as in the three-dimensional case 
(cf. to figure 4). The presence of the spin variables causes the critical value of k,2 to be 
shifted towards smaller values and the shift is largest when the spin system is critical. 
The nature of the transition between the two phases of the geometrical system seems 
to be the same at all points on the line k,2 = k-2{/3). There is a small increase of the 
height of the peak of the correlator (i?^) — (i?)^ as /3 — )• /3c but otherwise the data is 
compatible with the only effect of the matter field being a shift in k,2. In particular no 
significant change of the value of {R)[k,2) away from its pure gravity value is seen for 
any value of /3. 

Our partition function when we couple Ug Gaussian fields, . . . (f>^^ ^ to our geo- 
metrical system reads 

(4.239) 

where i and j label the 4-simplexes of our manifolds and where the (5-functions are 
introduced in order to eliminate the translational mode. We note that coupling a 
Gaussian field to simplicial gravity does not require the introduction of an additional 
coupling constant. A coupling constant in front of the Gaussian action can be absorbed 
in a redefinition of the cosmological constant by an appropriate rescaling of the scalar 
field. This of course refiects the fact that a Gaussian scalar field is automatically 
critical in the infinite volume limit. As in the case of the Ising spin system we see no 
drastic effects when Gaussian fields are coupled to the geometry. The critical value of 
is again shifted towards smaller values but the nature of the phase transition of the 
geometrical system seems unaltered. The height of the peak of the correlator (i?^) — 
{Ry is somewhat larger than for pure gravity, though. There are some indications of 
an enhancement of the interaction between matter and space time when the number of 
matter fields is increased. We have considered up till four Gaussian fields. The larger 
the number of Gaussian fields, the larger is the shift in ^2 the more pronounced 
is the peak in (i?^) — (i?)^- However the behaviour of (i?) as we pass the transition in 
geometry seems unaltered. In particular the value of (i?) at the critical points of the 
coupled systems is not significantly different from its value at the critical point of pure 
gravity. In the connection with the change in geometry for our space time manifolds 
we see a change in the expectation value ((^")^). The value of ((^")^) as well as the 
fiuctuations of this quantity increase when we enter the elongated phase. As in the 
3-dimensional case a detailed investigation of a possible change in critical indices for 
matter fields or geometry requires more extensive numerical studies. 
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4.4.3 Higher Derivative Regularization 



So far the space time part of the action for our models of dynamicaUy triangulated 
gravity has been given by the pure Einstein Hilbert action. In this section we will 
briefly describe the effects of adding to the Einstein Hilbert action a term. In the 
continuum formulation the introduction of a term (with positive coupling constant) 
would lead to a regularization of the path integral. In dynamically triangulated gravity 
the path integral involving only the Einstein Hilbert term is already regularized (for 
fixed space time volume). The purpose of introducing a R^ term in the dynamically 
triangulated model is to test whether this extension of the coupling constant space 
allows us to reach other critical points which might prove especially well suited for 
defining a continuum theory of quantum gravity. What we have in mind is of course a 
unmistakable second order phase transition point where the curvature scales to zero. 

To define simplicial analogues of quantities like fj^^ d'^^^R^, p > 1, one can ob- 
viously not maintain the original idea of Regge that the curvature is (5-function like 
having support only on the ( J — 2)-simplexes. A prescription for modifying the philoso- 
phy of Regge so that expressions like d'^^^R^ , p > 1 make sense has been given by 
Hamber and Williams [10]. The idea consists in viewing the curvature Rna-2 normally 
associated with a given [d — 2)-simplex, ^^-2, as being distributed in a J-dimensional 
volume around nd-2- One assigns a volume element Vd{nd-2) and a curvature density 
Rna-2 /^d{n-d-2) to each [d — 2)-simplex, ^^-2, in a J-dimensional simplicial manifold 
and calculates integrals involving powers of the curvature according to 

JM [Vd{nd-2)} 

Eor a J- dimensional simplicial manifold consisting of regular simplexes one has Rna-2 — 
(cd — o[nd-2)) and Ki^_2 = o[nd-2) / Kd+1,2 where we note that Kd+1,2 is the number of 
( J — 2)-dimensional sub-simplexes in a J-dimensional simplex. With these assignments 
and the prescription (4.240) for integrating over the space time manifold one gets 

/ d'CVaR' E oind-2)\"^^^^"^y Kd+1,2. (4.241) 
Jm [ o{nd-2) J 

Eurthermore the relations (2.7) and (2.10) are maintained. In the continuum formula- 
tion the addition of a R^ term to the Einstein Hilbert action results in the favourization 
of space time manifolds with zero local curvature. This statement does not immediately 
translate to dynamically triangulated gravity. In dynamically triangulated gravity it 
is not possible to have zero local curvature for d > 2. The addition of a R^ term like 
the one appearing on the right hand side of (4.241) does suppress manifolds with large 
numerical values of the curvature. However, the discretized R^ term takes its minimum 
value at i? = —0.46 not at i? = (where R was defined in equation (4.230)). Studying 
dynamically triangulated gravity with the lattice R^ term should not be considered 
as studying Einstein Hilbert gravity with exactly the term d'^^^R^ added but as 
exploring a certain class of models with higher derivative terms. Let us briefiy review 
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what has been learnt from such studies. The details can be found in reference [123]. 
In this work we consider 4- dimensional simplicial gravity with the following action 



fixed the value of h and looking for critical behaviour when varying k,2. The maximal 
value of h considered was h = 20. It is in practise impossible to increase h further 
because the acceptance rates in our Monte Carlo updatings become too small. For all 
values of h considered we encounter at a certain value ^2 of '^2 a transition between 
a hot and a cold phase of the geometrical system and in all cases the transition in 
geometry is accompanied by a peak in the susceptibility x('^2)/^4 the height of which 
tends to grow somewhat with volume. For not too large values of h the nature of the 
transition can not be distinguished from the nature of the transition encountered for 
h = 0. The value of (R) at the critical point does decrease with h but the change 
is very small. For h > 12 the peak in the correlator (i?^) — (i?)^ disappears. This 
might indicate that some new kind of transition takes place. However, (i?)(/€2, h) still 
decreases very slowly with h. Actually the data does not even allow us to conclude 
that (i?)(/€2,^) — > as /i — )• oo. Hence adding a higher derivative term to the action 
does not cure the problems encountered for pure gravity. 

Four-dimensional simplicial gravity with a term has also been studied in the 
static triangulation approach [122]. Here one finds that the first order phase transition 
encountered for pure gravity changes into a second order one when the value of h gets 
sufficiently large. Furthermore the value of (i?) at the critical point is always zero. This 
situation is very different from the one encountered for the dynamical triangulations. 
As in the 3-dimensional case the reason for this difference remains unclarified. 





The restricted coupling constant space of the model K4 = k,1[k,2^ h) is scanned keeping 
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